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Abstract: The triple differential spectrum in B — > X^ly is computed by Dressed Gluon 
Exponentiation (DGE). In this framework the on-shell calculation, converted into hadronic 
variables, can be directly used as an approximation to the meson decay spectrum, without 
involving a leading-power non-perturbative function. Sudakov resummation for the fully 
differential B — > X^lv width is formulated in moment space, where moments are defined 
using the ratio between the lightcone momentum components of the partonic jet /p~ 
and the hard scale is p~ . In these variables the correspondence with the B — > Xs'^ case 
is transparent. The Sudakov exponent is known to next -to-next-to-leading logarithmic 
accuracy. Further constraints are put on its Borel sum using the cancellation of the leading 
renormalon ambiguity and the absence of the next-to-leading one, which was proven in 
the large-/3o limit and assumed here to be general. Based on the resummed spectrum, 
matched to the fully differential NLO result, we calculate the event fraction associated 
with experimental cuts on the hadronic mass (or the small lightcone component) as well 
as on the lepton energy. Finally, we extract from recent measurements by Belle and 
analyze the theoretical uncertainty. 
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1. Introduction 

A central issue on the agenda of the B factories is the precise determination of the CKM 
matrix element V^b- Without undermining the significance of exclusive hadronic final- 
state analysis, measurements of the inclusive branching fraction of charmless semileptonic 
decays, B — > Xylv, provide the safest and most accurate determination of \Vuh\ 
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Theoretically, because of the inclusive nature of the measurement, the calculation of 
the width F (5 — > X^W) in units of \Vubf relies primarily on QCD perturbation theory. 
The final state, which in reality is composed of light hadrons, can be safely described in 
the course of the calculation by light quarks and gluons. Moreover, the dependence of the 
width on the internal dynamics of the decaying B meson is power suppressed. The Operator 
Product Expansion (OPE) implies that the total decay width of a B meson is equal to the 
total decay width of a hypothetical on-shell b quark up to ©(A^/m^) corrections, which 
can be expressed as forward matrix elements of certain local operators whose numerical 
values are determined based on other measurements. These corrections are very small. 

The main obstacle in extracting \Vui,\ from inclusive semileptonic decays is the need 
to eliminate the huge background due to B decay into charm. The most effective way to 
maximize the rate while discriminating the charm background is to select events where the 
hadronic system has small mass, Mx ^ 1-7 GeV, sec e.g. [5,7-9]. 

Because of the cuts, detailed theoretical understanding of the differential distribution is 
essential. In the small-Mx region the hadronic system is most likely jet-like. Working with 
lightcone coordinates, P''" < P" < Mb with M|- = P'^P~, the small Mx region directly 
corresponds^ to the region where the smaller lighcone component gets small, P+ — > 0, so 
it involves highly non-local correlation functions with large lightcone separation, y" — > oo, 
where y~ is the Fourier conjugate to P~^ . At leading power in 0{A/mh) the differential 
width involves the full lightcone momentum distribution function of the b quark inside the 
meson, the so-called "shape function" [10,11]. On the non-perturbative level, not much is 
known about this distribution: it is difficult to compute it on the lattice and although in 
principle it can be measured using the photon spectrum in rare radiative decays B — ^ Xg^, 
realistic measurements constrain just the first few moments of this function [12,13]. 

As it stands, the theoretical uncertainty in extracting \ Vub\ from inclusive semileptonic 
decays is dominated by the uncertainty in estimating the effect of cuts. A great variety of 
experimental cuts has been proposed, which either aim at reducing the sensitivity to the 
unknown momentum distribution of the b quark in the meson by cutting out the phase- 
space regions where it is important [14,15], or alternatively, aim at making the relation 
with the B — > Xsj spectrum as direct as possible [16-19]. While the comparison between 
results for \Vub\ obtained using different cuts is useful for controlling experimental as well 
as theoretical uncertainties, better measurements can be achieved if cuts would be chosen 
based on purely experimental considerations. Precise theoretical predictions for the fully 
differential width are therefore urgently needed. 

1.1 The applicability of perturbation theory 

In perturbation theory, the fully differential width is currently known [20] to the next- 
to-leading order (NLO), which includes tree level as well as one-loop, i.e. 0{as), virtual 
corrections, where the hadronic system is represented by a single on-shcll u quark, as 
well as bremsstrahlung contributions where the hadronic system is composed of a u-quark 
and a gluon, and therefore has a non-vanishing mass. It has long been recognized that 

^The region where both lightcone components are small is power suppressed. 
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fixed-order calculations do not provide even a qualitative description of the differential 
width in the small-Mjf region, where multiple emission of soft and collinear gluons is 
important [21-31]. In terms of partonic^ lightcone coordinates, p'j' < pj < nib, the small 
Mx region is characterized by a large hierarchy, pj/pj 1. Higher-order corrections 
containing Sudakov logarithms of the ratio /pj are large, and they must be resummed 
to all orders to recover the characteristic Sudakov peak at small Mx- 

According to the factorization properties of the fully differential width in inclusive 
decays [21-31], Sudakov logarithms exponentiate in moment space. We shall work here 
with moments n defined^ with respect to powers of 1 —Pj'/pJ, as in Eq. (3.18) below. The 
region of interest, of small p^, is probed by high moments, n — > oo. The logarithms. 
Inn, which become large in this limit, originate in two distinct subproccsscs of different 
characteristic scales, which factorizc in this moment space. The first is the final-state jet 
having a constrained invariant mass squared of 0{m'^/n) [30-32], which coincides with 
the Sudakov factor of deep inelastic structure functions. The second corresponds to the 
quark distribution in an on-shell heavy quark, the perturbative analogue of the quark 
distribution in the meson [30,31,33], and describes soft radiation, \k\ ~ 0{mb/n), from the 
nearly on-shell heavy quark prior to the decay. 

A major stumbling block, which essentially prohibits the straightforward application 
of Sudakov rcsummation to inclusive B decay phenomenology, is that the scales involved 
are almost too low for perturbation theory to apply. While the b quark mass is heavy 
compared to the QCD scale A, in the Sudakov limit, n — > oo, the jet mass scale is much 
lower than and the soft scale |A;| ~ 0{mi,/n) characterizing the quark distribution is yet 
much lower. Already at moderately high moments, n ~ 10, the latter is dangerously close 
to A where the physics in non-perturbative. On these grounds it is often argued [34,35] that 
Sudakov resummation is useless in this case. Conversely, when Sudakov resummation is 
applied it is usually supplemented by an external infrared cutoff along with parametrization 
of the non-perturbative quark distribution function below this scale [16, 17]. 

While there are many ways to partially bypass the need to resum corrections on the 
scale \k\ ~ 0{mh/n), we argue that this resummation is in fact largely under control and 
that it is very useful. Ref. [31] and the present investigation show that definite predictions 
emerge when the available information on the large-order behavior of the perturbative 
expansion is taken into account. The sensitivity to low momentum scales, and specifically 
to the meson structure, turns out to be smaller than it superficially appears to be, and 
consequently the predictive power of perturbation theory is higher than what one might a 
priori expect. 

Infrared sensitivity appears in the moment-space Sudakov exponent through infrared 
renormalons [32, 33, 36-41]. These generate divergences and invalidate the logarithmic 
accuracy criterion [30,31], but they definitely do not prohibit using Sudakov resumma- 
tion altogether: the Sudakov exponent needs to be regarded as an asymptotic series, and 
summed up to all orders by regularizing the renormalons in a systematic way. This is the 

^The relation with the hadronic Ughtcone coordinates is pf = —A where A = MB — rrib, see Eq. (4.1). 
^There are other viable definitions for moments; an alternative is considered in Appendix A, see 
Eq. (A. 12) there. 
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essence of Dressed Gluon Exponentiation (DGE). The application of this approach to inclu- 
sive decays is particularly advantageous since the leading infrared renormalon ambiguity, 
which dictates the divergence of the series in the Sudakov exponent, cancels exactly with 
the pole-mass ambiguity that is carried by kinematic power corrections associated with 
the conversion from partonic to hadronic variables. In this way, the on-shell calculation 
becomes directly useful for phenomenology, in spite of the divergence of the series (or the 
ambiguity) that appeared at the intermediate stage. 

In the present paper we provide precise theoretical predictions for the fully differential 
width in semileptonic decays using DGE. Our calculation makes maximal use of Sudakov 
resummation and the underlying renormalon structure. In this way we effectively minimize 
the role of the unknown non-perturbative component of the quark-distribution function 
and the associated uncertainties. 

The method of DGE was already proven successful in the application to the photon 
energy spectrum in B — > Xsj decays [31,41]. Specifically, the first two central moments 
defined over the range > Eq, were computed in [31] in an essentially perturbative fashion 
(see below) without involving any parametrization of the leading-power quark distribution 
function. These predictions were later compared [41] with experimental measurements from 
Babar [42] over a wide range of cuts, Eq = 1.9 to 2.26 GeV, finding very good agreement. 

1.2 Comparison with alternative approaches 

It is worthwhile re-examining from the present perspective the conventional approach to 
inclusive spectra, which is based on parametrizing the leading-power non-perturbative 
"shape function" [10,11] in analogy with deep inelastic structure-function phenomenology. 
Let us summarize the qualitative differences. 

In Refs. [10, 11, 17] one begins by considering directly the kinematic region where the 
small lightcone component of the hadronic system is comparable to the momentum of the 
light degrees of freedom in the meson, 0{A). In the — > oo limit the long-distance 
interaction is captured by a single lightcone distribution, the "shape function" . No attempt 
is made to compute this function, which is assumed non-perturbative from the beginning. 
The "shape function" is parametrized and the parameters are fixed by a fit to the measured 
B — y Xg'y spectrum (see e.g. Sec. 4 in Ref. [17]). Owing to the universality of this function, 
the result can be directly used in semi-leptonic decays. There is of course some bias owing 
to the functional form assumed. In practice this can be dealt with by varying the function, 
or better, by relating directly the semileptonic distribution to the B — > Xg^ data by 
weighted integrals [18,19]. 

Our approach is more ambitious: we wish to compute the spectrum. Relying on the 
fact that the heavy quark in the meson is rarely far off shell, we compute the on-shell decay 
spectrum and use it as a first approximation to the meson decay spectrum. In this approx- 
imation the quark distribution in the meson is replaced by its perturbative analogue, the 
quark distribution in an on-shell heavy quark — see Ref. [33] for the precise definition. A 
key element in this approach is the systematic treatment of infrared renormalons. Most 
importantly, this includes a systematic definition of the on-shell heavy quark state at the 
level of resummed perturbation theory. Having full control of the large-order behavior 
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of the expansion, such definition is provided by Principal Value Borel summation. Upon 
consistently using this definition in both the Sudakov exponent and the quark pole mass 
(which enters the computed spectrum through the conversion to hadronic variables) any 
0{A/mh) ambiguities are avoided. In this way resummed perturbation theory yields def- 
inite predictions for the on-shell decay spectrum. The latter provides an approximation 
to the meson decay spectrum that does not require any 0{A/mb) non-perturbative power 
corrections. Having established that, we study higher power corrections that constitute the 
ratio between the moments of the physical spectrum and the computed one, which originate 
in the dynamical structure of the meson. While these corrections gradually increase with 
the moment index n they have but a small effect on the global properties of the spectrum. 

It should be emphasized that in reality the spectrum in the immediate vicinity of the 
endpoint, P"*" ~ 0, cannot be accurately described by any of these approaches for two 
reasons: first, the number of relevant "shape function" parameters increases, and second, 
the transition to the exclusive region, namely the hadronic structure of the final state, is 
not addressed. 

Recent formulations of the "shape function" approach, e.g. [17], employ the tools of 
factorization [21, 22] and can therefore make some use of Sudakov resummation. Nev- 
ertheless the resummation there effectively concerns the jet function only, as the "shape 
function" is defined at an intermediate scale of order of the jet mass. Moreover, the need 
to separate "jet" logarithms from "soft" ones requires a complicated infrared cutoff to be 
introduced at any order in perturbation theory. The final answer is affected by residual 
scale and scheme dependence. 

In contrast, in our formulation Sudakov resummation is applied also to the quark- 
distribution function. We thus make full use of the available perturbative information. 
Moreover, the on-shell decay spectrum we compute is itself renormalization group invariant: 
it is absolutely free of any factorization scale or scheme dependence; at the level of the 
Sudakov factor, it is also free of renormalization-scale dependence owing to the additional 
resummation of running-coupling effects we perform. 

The conclusion from Ref. [31] and the successful comparison with experimental data [41] 
is that when using DGE, the on-shell calculation, b — > Xg'j, converted into hadronic vari- 
ables, directly provides a good approximation of the B meson decay spectrum. Notably, 
in contrast with any fixed-order approximation, it has approximately correct physical sup- 
port properties at large E^. This implies, in particular, that in this framework the non- 
perturbative component of the lightcone momentum distribution function is small, and at 
the present accuracy can even be neglected. This conclusion may be surprising: as dis- 
cussed above, in other approaches the paramctrization of this non-perturbative distribution 
is absolutely essential. In order to understand the difference one should recall that 

(1) The separation between the perturbative and non-perturbative components of the 
lightcone momentum distribution function depends on the convention adopted. In 
our approach, the calculated perturbative component is nothing but the quark distri- 
bution in an on-shell heavy quark. It requires renormalization owing to logarithmic 
ultraviolet singularities, however, no infrared cutoff is needed. This distribution is 
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infrared and coUinear safe and therefore well-defined to all orders in perturbation 
theory. It requires of course Sudakov resummation, owing to the large hierarchy 
between the hard scale and the soft scale rrib/n. In contrast, at the power level 
the Sudakov exponent is infrared sensitive: its Borel sum has ambiguities scaling 
as powers of (nA/mfe). In our approach these are regularized using the Principal 
Value prescription — this is how we define the perturbative component of the quark 
distribution function. 

(2) Beyond the currently available logarithmic accuracy at which the quark distribu- 
tion is known (next-to next -to-leading logarithmic accuracy [33]), the DGE calcu- 
lation makes use of additional information, which is important in constraining the 
corresponding Borel sum. Most importantly, the exact cancellation of the leading 
{u = 1/2) renormalon ambiguity with kinematic power corrections involving the pole 
mass [30] is used to fix the corresponding residue [31]. In addition, the absence of 
the next-to-leading renormalon (u = 1), which was proven in the large-/3o limit [30], 
is assumed here to be general. Thus, the dependence on the (Principal Value) pre- 
scription is eventually restricted to higher renormalons {u > 3/2), which have limited 
influence on the spectrum — see Ref. [31] and below. 

Following the results of Refs. [31,41], our working assumption in this paper is that the 
properly-defined quark distribution in an on-shell heavy quark provides a good approxi- 
mation to the quark distribution in the meson. In reality the two differ of course, and it 
is important to quantify the power corrections that distinguish between them, a task that 
will require theoretical as well as experimental input. This, however, is beyond the scope 
of the present paper. 

1.3 The task and the strategy 

The main task on which we embark here is to provide a reliable calculation of the partial 

decay width with experimentally-driven cuts. Priority is given to cuts that maximize the 
rate, in particular [7] the charm-discriminating cut on the hadronic mass P^P^ < = 
(1.7 GcV)^ with an additional, experimentally unavoidable, mild cut on the charged Icpton 
energy Ei > 1 GeV. Our general strategy to extract \Vub\ from data is based on separately 
computing: 

• The total charmless semileptonic width Ftotai {B — > Xulv) in units of using 
the available NNLO result [43] and renormalon resummation. 

• The effect of experimental cuts as the event fraction 

^ _ T{B — > Xylv restricted phase space) ^, 
V{B — > Xylv entire phase space) ' 

which we compute using DGE and match to the fully differential NLO result in 
moment space. 
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The partial branching ratio measured by the B factories can then be compared to 

AB(S — > Xuli> restricted phase space) = rsFtotai (-B — ^ X^liy) -Rcut- (1-2) 

The advantages in spHtting the calculation of the partial width this way are: 

• The event fraction Rent is free of the 0{A/mh) renormalon ambiguity associated 
with the overall factor where rrih is the quark pole mass. Note that 0(nA/mb) 
cut-related renormalon ambiguities are present in the perturbative result for i^cut, 
but their effect is restricted to the moment-space Sudakov exponent. Therefore, the 
perturbative expansion of the hard matching coefficient for .Rcut is renormalon free 
and truly dominated by hard scales at higher orders. 

• The total charmless semileptonic width can be computed with NNLO accuracy and 
in a way that renormalon-related 0{A/mh) power effects explicitly cancel. One then 
obtains the total width with theoretical uncertainty as low as ~ 6%. Such accu- 
racy cannot be achieved in a direct calculation of the partial width in the restricted 
phase-space since (1) the differential width is known in full only to NLO accuracy; 
and (2) while Sudakov resummation can be used to take into account parametrically- 
enhanced corrections associated with the phase-space restriction, some residual un- 
certainty of both perturbative and non-perturbative nature remains. Quantitative 
estimates are provided in Sec. 4.2. 

The paper is divided into three main sections. In Sec. 2 we compute the total semilep- 
tonic width using the NNLO result [43] and renormalon resummation. In Sec. 3 wc present 
the results for resummation of the triple-differential width and work out NLO matching 
formulae in moment space. Finally, in Sec. 4 we use these results to compute the differ- 
ential and partially-integrated width in hadronic variables, taking into account the mass 
difference between the B meson and the b quark. In this section we study the theoretical 
uncertainty in i?cut and extract \Vuh\ from recent measurements by Belle [7]. 



2. The total chcirmless semileptonic decay width 

According to the OPE the total B — > Xylv decay width is given by 



r,„,, [B XJu) = ^£i^^F(a,(m,)) (2.1) 



1927r3 
G%\Vubfml 
1927r3 



I 1 1 

TT \ TT / 



ml 



where the perturbative expansion has been computed to NNLO [43] and the matrix ele- 
ments controlling the O (A^/m^) terms are not large and amount to ~ 1% correction. 

It is well known [44-46] that while the decay width receives just small 0(A^/m^) 
power corrections, the perturbative expansions of (a) the ratio between the pole mass 
and short-distance masses (e.g. m^'^^); and (b) the function F{as{mi,)) in Eq. (2.1) have 
large 0(A/mh) renormalon ambiguities owing to infrared sensitivity. Furthermore, the 
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first few terms in these expansions show poor convergence; the asymptotic behavior sets 
in early. A practical problem one needs to address is how to optimally use the known 
coefficients in these expansions to get a reliable estimate of Ttotai in Eq. (2.1). Often 
(e.g. [17,34]) this problem is solved by resorting to a specific mass scheme where the mass 
and the perturbative expansion of the width are separately renormalon free, and hopefully 
converge well. Our approach is different. We use the pole mass itself and regularize the 
leading renormalon using Borel summation^, explicitly using the cancellation in Eq. (2.1). 

In Ref. [31] we studied the mass ratio and expressed it as a Borel integral where the 
Borel function is written as a bi-local expansion [48-50] 



nib 



= 1 + 



' f'OO "-max 

'0 ^0 



+ 



Q 

(1 - 2z) 



1 + ^0^(1-2^)^= 



k=l 



(2.2) 



where 5 = Pi/Pq, cj. depend on the coefficients of the P function (see Appendix B in [31]), 
and the coefficients bi in the regular part of the expansion can be determined through order 
62 knowing the N^LO result [51] for the mass ratio. Here both the exact analytic structure of 
the singularity [52] and the value of the residue q fixing the large-order asymptotic behavior 
of the expansion arc used. The latter can be determined with ^ 3% accuracy [31,50,53,54] 
using the N'^LO result for the mass ratio. Hence, given a definite integration prescription 
in the complex z plain to avoid the z = ^ singularity, Eq. (2.2) facilitates an accurate 
calculation of both the 0(1) real and the 0{A/mb) imaginary parts of the mass ratio. 
With the central values of the residue from Ref. [31], 



CpqiNf = 3)/7r = 0.560; CpqiNj = 4)/7r = 0.536, 



one obtains: 



nib 



mt 



= 1.161 - 0.0934 z 



mb 



Nf 



ml 



= 1.164 - 0.0768 z. 



(2.3) 



(2.4) 



Nf 



respectively. The variation in the number of light fermions is used to estimate charm mass 
effects. Here we chose to define the Borel sum in Eq. (2.2) with the integration contour 
going below the real axis avoiding the singularity at z = ^. This fixes the imaginary part 
of the mass ratio (the physical real part is prescription independent). We evaluated the 
integrals using: 

-^/^ A „ c-'r ZQ 

r(i + y) Jo 



f 

Jo 



dz- 



(-4) 



dl 



TT 



1 _ M Til + v) suii-Kv 



e-^°l^{-ZQlAY 



+A 



i^e 







u 



1 



(2.5) 



*The possibility to determine the normalization of leading renormalon residues using the structure of the 
singularity and the first few orders in the perturbative expansion was considered by T. Lee (in a different 
context) already in [47]. It was further observed that this information can be used to systematically 
construct a bi-local expansion for the Borel transform [49] . Then Pineda found [50] that in the case of the 
pole mass the residue can be accurately determined, and used it to subtract the corresponding divergence. 
More recently, it was shown that the Principal Value Borel-resummed pole mass can be directly computed 
from the bi-local expansion [31,48] and used as an alternative to short-distance or renormalon-subtracted 
mass definitions. 
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where we assumed that zq = ^ + ie where e is an infinitesimally small positive number, and 
Mn^u{a) is the Whittaker M function, 

1 



+ v — H 



\ + 2v 



,a 



The next crucial observation is that the cancellation of the leading renormalon ambi- 
guity in Eq. (2.1) implies that 



lm{F^/'\as{mb))} 
Re{FV5(a,(m6))} ~ 



Im ' 






Re< 


nib I 





0(AM), 



(2.6) 



and that the observable itself can be computed with 0(A^/m^) accuracy using the real 
parts (or Principal Value) of the mass ratio and F^/^: 

rrn,Fy'iasim,)) = x Re x Re {FV5(a^(^,))} . (2.7) 

Eq. (2.6) implies that the ambiguity of F^/^{as{'mb)), just like that of the mass ratio [52], 
is a pure power term (i.e. a number times A/ mi,), not modified by logarithms. This is an 
exact result. 

The structure of the Borel singularity is therefore: 



/-t poo Jmax 
PO Jo ~~r. 



bi + 



(2.8) 



(1-2^)1+^5 

where the state-of-the-art knowledge of the expansion of F{Q.s{m}))) in Eq. (2.1) (NNLO 
accuracy) allows determination of ho and hi (so jmax = !)• Evaluating the Borel sum in 
Eq. (2.8) and using Eq. (2.6) and Eq. (2.4) we obtain 

g|jv,=3 = -1-191 ; ^liv,=4 = -1-146 , (2.9) 

respectively. This yields: 

F^/^{as{mi,)) = 0.933 + 0.0750 F^/^{as{mi)) = 0.928 + 0.0612 z , 

(2.10) 

(2.11) 



iVf=3 



N(=A 



respectively. Using Eq. (2.10) and Eq. (2.4) with 

= mf^{mf^) = 4.19 ± 0.05 GeV 
we find 

mfF(a,(mb))|^^^3 = 1926 ± 117 GeV^ mfF(a,(mfe))|^^^4 = 1901±116GeV^ 

(2.12) 

Thus, in units of the total semileptonic decay width is 

^^r,„,,, {B XJv) = 66.5 ± 4ps-\ (2.13) 



IK 



ub\ 



where the error is dominated by the uncertainty in the value of the short-distance mass 
in Eq. (2.11). The central value as well as the error estimate are consistent with previous 
studies, see e.g. Ref. [55]. 
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3. Resummation of the triple differential width 

3.1 Partonic kinematics and triple differential width at NLO 

In order to describe the triple differential width in b{pb) — > Xu{pj) l{ki) v{ki,) we define 
the following kinematic variables: the charged lepton energy fraction, 

X = 2pb- k/ml, (3.1) 

the lightcone variables defined by the partonic jet 

p = pprrib = {Ej - \pj\)/mb, 

A = p-/mb = {Ej + \pj\)/mb, (3.2) 

and their ratio 

P_Pl_E^-jp^ 

x-pj- E, + \p,y ^'-'^ 

which is the exponent of the rapidity. As we shall see below the structure of the Sudakov 
exponent and the matching is particularly simple when expressed in moments of 1 — r. 
Using the variables (A, r, x) the phase-space integration is: 

The perturbative expansion of the triple differential width takes the form: 

1 dr(A,r, X 



To dXdrdx 



= V{\x)5{r) + R{\r,x) (3.5) 



where 



and 



= 192.3 " (^■^) 
F(A,x) = ^o(A,x) + ^^^^^«;i(A,a;) + ---, (3.7) 

TT 

i?(A,r,x) = ^^^^^MA,r,x) + ..., 

TT 

correspond to virtual and real corrections, respectively. 

The 0[a.^ coefficients (NLO) are known since long [20], however, they were computed 
there in terms of the hadronic mass and energy variables which are less suited for resum- 
mation as compared to the lightcone variables introduced above. Using (A, r, x) the virtual 
coefficients are: 

u;o = 12 (2 - X - A) (x + A - 1) (3.8) 

5 vr^ 



u;i = 12 (2 - X - A) (x + A - 1) 

+6(2A + 2x-5) ln(A) (x + A - 1 



----Li2(l-A) 
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and the real-emission coefHcient is 
61nr 



bin T ( \ 

ki{X,r,x) = — —r.{x-2 + X + rX)(rX{l-X) + x + X-l] x 

r [r — ij^ \ / 

2 + 2r^ + r^X'^-2r^X-6r'^X + r'^X'^ + 2r'^ + 2r-2rX 
3 



(3.9) 



+ 



(A + rA-rA^ + x-l)(a;-2 + A + rA) 



+ 



+ 



r(l — r) 

Ar^ X'^ - lOr'^ X + 7 r'^ + 2r - lOr X + 7 

6 (-2 r2 A + 2 r2 + A2 - 2 r A + 4r + 2) (-1 + r A) (-1 + A) A ln(r) 

(l-r)2 

3 (r^ A^ - 8 r A + 8 r + r A^ + 8) (-1 + r A) (-1 + A) A 

1 — r 

6(x — l + rA)(A + x — l)lnr 



(l-rY 



-14r + 42r^A-2rA^-2r^A^ + 16r^A 



+ 2r* A - 22r2 A^ + 16r A - 4r'' + 2 A - 22r'' A^ - 14r^ + A'^r'' - 4 



3 \2 



+ Ar^X^ + r^ X^ 



3(x-l + rA)(-l + A + x) 
(l-r)3 



Gr^ A^ + llr^A- IGr^ A^ 



+ 67 A - 64r^ A^ - 22r^ + 6r^ A^ - 16r A^ - 28r + 67r A - 22 + 11 A 



The r — ^ singular (non integrable) terms in ki, namely 



fep(A,r,x) = -12(2-a;-A)(x + A-l) 



Inr 
r 



4 \r 



-Wo 



Inr 

r 



-(- 

4 Vr 



are regularized as ()* distributions, i.e. 

'' ° drF{r) (^^^ =F(0)lnro + J " dr(^F{r) - F{0)^^, 







/•'■o /lnr\ 1 /'''° / 

J drF{r)i^—j = F{0)-ln^ ro + J dr(F{r) - F{0) 

where F{r) is a smooth test function. The remaining, regular part of ki 

(A, r, x) = ki (A, r, x) - (A, r, x) 



Inr 
r 



(3.10) 

(3.11) 
(3.12) 

(3.13) 



is integrable, so it does not require any regularization. 

As usual, at 0{as) and beyond the separation into real and virtual terms is not unique. 
In order to compare the result quoted here with NLO expressions for the triple differential 
width in terms of other kinematic variables (e.g. in [20]) one must take account of the 
fact that terms proportional to 6{r) are contained both in the ()* distributions and in the 
virtual terms and, depending on the variables chosen, they may be split differently between 
the two. 
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3.2 The Sudakov exponent 

Owing to multiple soft and collinear gluon emission, Sudakov logarithms, such as the 
r — > singular terms of Eq. (3.10), appear at any order in perturbation theory. Such 
terms spoil the convergence of the perturbative expansion and must be resummed. Sudakov 
logarithms in inclusive decay spectra are associated with two independent subprocesses [21, 
30,31,40] (see also Refs. [22-29]): the soft function, which is the Sudakov factor of the 
heavy quark distribution function [33] and the jet function, summing up radiation that is 
associated with an unresolved final-state quark jet of a given invariant mass. This function 
is directly related to the large-x limit of deep inelastic structure functions; see Ref. [31,32] 
for further details. 

Singular terms in the large— /?o limit 

Ref. [30] has generalized the concept of Sudakov resummation in inclusive decay spectra 
beyond the perturbative (logarithmic) level. It has been shown that, when considered 
to all orders, the moments corresponding to each of the above subprocesses contain in- 
frared renormalons and therefore certain power corrections exponentiate together with the 
logarithms. 

The calculation in Ref. [30] was based on evaluating the real-emission diagrams with 
a single dressed gluon using the Borel technique: the dimension of the gluon propagator is 
modified: l/{-k^) — > l/(-/c2)^+", and in this way one gains control of the physical scale 
for the running coupling. Ref. [30] used the leptonic and partonic invariant mass variables 
z = q^/ml and 1 — ^ = pj/mf, that are related to the lightcone variables of Eq. (3.2) by: 

z = {1 - p){l - X); e = l-Ap, (3.14) 

and obtained the following result for the ^ — > 1 singular terms in the triple differential 
width in the large-/3o limit: 



1 dr{^,z,x) 



To dzdxd^ 
+ 



= l2{x- z){l + z-x) X {5{l-0+ (3.15) 

e— ^1 



Cf rdu fA^Y 



Bs{u){l - zf^ (1 - e)-'-'" - Bjiu) (1 - 0' 



+ 



where the dots stand for terms that are subleading in Po or are suppressed by powers of 
(1 — ^), and the Borel functions corresponding to the two Sudakov anomalous dimensions 
are: 

Bs{u) =eh(l-u) + O(l//3o), 

Bjin) = \ et" {-^ + -i-) ^ + Oil/M. (3.16) 

2 \1 — U 1 — U/2 J TTU 

Here and below we use the scheme-invariant Borel transform [57] where T{u) is the Laplace 
conjugate of the coupling. In the large-/3o limit (one-loop running) T{u) = 1. Below we 
work in the full theory and define T{u) as the Laplace conjugate of the 't Hooft coupling, 
see Eq. (2.18) in Ref. [31]. 
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Eq. (3.15) reveals the physical scales that control the running coupling within the soft 
and the jet subprocess, respectively: the soft scale is \k\ ~ mb{l —$,)/{! — z) while the 
invariant mass squared of the jet is A;^ ~ m^(l — ^). 

Exponentiation 

Sudakov rcsummation can be done using different variables. One natural possibility is 
to resum Sudakov logarithms of 1 — corresponding to the limit where the hadronic 
invariant mass gets small, for fixed leptonic invariant mass z. This avenue is followed in 
Appendix A. Here instead we choose to revert to the lightcone variables of Eq. (3.2), where 
the exponent and the matching procedure are simpler. The comparison with Appendix A 
provides a useful consistency check. 

The Sudakov limit — > 1 corresponds to the phase-space limit where the smaller 
lightcone momentum component p tends to zero. The large component A can either be 
small or of 0{1); fortunately, the region where also A is small is power suppressed in the 
Born-level weight (see Eq. (C.2)) and is therefore unimportant. Consequently, Sudakov 
resummation can be readily applied to resum logarithms of p or of r. The latter possibility 
is particularly attractive. The singular terms for r — > in the large-/3o limit take the 
form: 

1 dT{X,r,x) 



dXdrdx 



= Wo { S{r) + 





so both the soft and jet functions depend just on r, except for an overall dependence on A 
that scales the argument of the coupling. It is this dependence that makes the difference 
with the radiative decay case [31], which simply corresponds to the substitution of A = 1 
in Eq. (3.17). 

Multiple emission of soft and collinear gluons is taken into account to all orders by 
exponentiation of the fully differential width in moment space. It is necessary to consider 
the fully differential width, since only then the hard configuration, from which soft and 
collinear gluons are radiated, is fixed^. The integration over the variables A and x that 
control the hard configuration itself can only be performed after exponentiation has been 
carried out. Exponentiation takes place in moment space because the phase space of soft 
and collinear radiation factorizes there. Defining moments by 



~rJo ^ ' dXdrdx ^'^■^^> 

n-^dr{X,r,x) ^^^^^^^^ 



dXdx 



-f 

To in 



dr (1 — r 



dXdrdx 



4-0 



the singular (non-integrable) terms for r — > of the form appearing in Eq. (3.17) generate 
terms that contain powers of Inn, while integrable terms generate 0{l/n) terms, that are 

^This point has recently boon pointed out in Ref. [23]. We note that the variable u introduced in that 
paper is identical to our variable r. 
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neglected at this stage. These terms will eventually be taken into account (to a given 
order) by matching the resummed result to a fixed-order expression. The resummed, 
triple differetial distribution can then be computed to all orders by the following inverse 
Mellin transformation: 

1 dr{X,r,x) _ [^+'^ dn dTniKx) _ ^yn (3 



Fq dXdrdx Jc~ioo ^Tri dXdx 

Computing the moments in Eq. (3.18) with Eq. (3.17) and exponentiating the terms that 
diverge in the large-n limit we get the following DGE formula: 



dr„(A,x) 



dXdx 



= [ViX,x) + AR^iX,x)] xexpj^ r^Tiu)(^) (3.20) 

large n [ Po Jo U \mfX^ J 

[Bs{u)r{-2u) (n^" - 1) - Bj{u)r{-u) (n" - 1)] | + 0{l/n) 

where the square brackets that multiply the Sudakov factor contain the virtual terms 

V{X, x) of Eq. (3.5) as well as Ai?°° that is defined as the n — > 00 limit of the difference 
between the moments of the full real-emission contribution -R(A, r, x) and the terms that 
are included in the exponent, at any given order. Ai?°° will be determined in the next 
section at 0{as) in the process of matching Eq. (3.20) to the full NLO expression; see Eqs. 
(3.35) and (3.37). 

Eq. (3.20) summarizes the exact, all-order structure of the Sudakov factor. Comparing 
Eq. (3.20) to Eq. (3.17), which was derived in the large— /3o limit, one notes that terms that 
are subleading in (3q arise from several sources: (1) the exponentiation; (2) the function T[u) 
summarizing the dependence on the two-loop (5 function coefficient (5i ; (3) the anomalous- 
dimension functions Bsiu) and Bj{u) that receive ©(l/Z^o) contributions starting at 0{u). 
The Sudakov factor in B — > Xg'y [31] can be recovered by substituting A = 1 in Eq. (3.20). 

The perturbative expansions of Bs{u) and Bj{u) have recently been determined [31,33] 
to O(n^), i.e. the NNLO; explicit expressions appear in Sec. 2.2 in [31], see Eq. (3.27) 
below. In the following we shall compute the Borcl sum in Eq. (3.20) directly, using the 
Principal Value prescription. Before doing so, however, it is worthwhile looking at the 
conventional approach to Sudakov resummation, where a logarithmic-accuracy criterion is 
used. 

Resummation with a fixed logarithmic accuracy 

Let us write the exponent in Eq. (3.20) as an integral of the Sudakov anomalous dimensions 
over the range of scales: 

'^^^ = [V{X,x) + /^R°-{\x)\ xexpi r^^ ^^-^)""'-^ (3.21) 



'^A{as{^Ji^)) + B{as{mlX\)) - V{as{mlX\'')) 
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where the anomalous dimension functions have the following expansions in the MS scheme: 

i2\ 



oo 
n=l 






oo 
n=l 






n=l 







duT(u) 



2 \ 



A2 



Bb{u), 



(3.22) 



where the known^ coefficients ^i,2,3) -81,2 and Di^2 are detailed in Eqs. (2.5) through 
(2.7) in Ref. [3f] and the all-order relation with the Borel functions Bs{u) and Bj{u) of 
Eq. (3.20) above is 



Bs{u) = B_a{u) -uBv{u), 
Bjiu) = Bj^{u)-uBb{u). 

The moment-space Sudakov factor of Eq. (3.21) can be written as: 



(3.23) 



where 



Sud(mbA, n) = exp O Qkir) ( ^ J 

U=o \ / 



k-r 



(3.24) 



T = 



TT 



-/3o Inn. 



A fixed-logarithmic-accuracy approximation is obtained by a given truncation of the sum 
over k in the exponent. The first three coefficients gk{T), which sum up the logarithms to 
NNLL accuracy, are: 



50 (r) 
9i{r) 



= ((1 - r) In (1 - r) - i (1 - 2r) In (1 - 2r)) 



A 

A 
1 



(3.25) 



In (1 - r) + In (1 - 2t) ) + Bi In (1 - r) - - L>i In (1 - 2t) 



+ 



+ 



A2_ 



- In (1 - r) + - In (1 - 2r) 

iln(l-2r)-iln(l-2r)2 + iln(l-r)2 + ln(l-r) 



TT 



®The coefficient A3 was recently computed in the course of the three-loop calculation of the Altarelli- 

Parisi evolution kernel [58]. B2 has been known for some time [32,59] based on two-loop calculations in dccp- 
inelastic scattering. D2 was recently computed [33] directly from the renormalization of the corresponding 
Wilson-line operator [60] at two-loops. Independent confirmation of this result has been possible owing to 
an exact all-order correspondence [33] between the heavy quark distribution and fragmentation functions 
in the Sudakov limit and a recent two- loop calculation of the latter [61]. 
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+Bi 7B 1 



1 



1 



1-T 
1 1 



+ Di7E 



l-2r 



A27E + 1 



+ 



1 

Pi 



+ 



1 



1-T l-2r 

In (1 - r) _ ln(l - 2r) 
r - 1 2r - 1 

In (1-T 



1-T 



+ 



B2 + 

1 



T-1 2t-1 



2 4t-2 



£'2 



T-1 



T-1 



_1 1_ 

' 4'^8t-4 2t-2 



1 I 5i + 

^3 



1 1 lln-2T + l 
- H \ ^ 

2 4t- 2 2 2t-1 



Pi 



1 / /In (1-T 



3 3 1 ln(l - 2t) 3 

T-1 '2 (t-1) 4 ~ 2 2t- 1 4 (2t - 1) 



+ 



+ fln(l-T)--ln(l-2T)-- —+ ^ 

V ^ ^ 2 ^ ' \ 2t-2 



T-4 



Po' 
In (1 - 2t) t 

2t- 1 



1 In^ (1 - t) _ In (1 - t) T 

2 T-1 T-1 
1 1 



1 _^ 1 In^ (1 - 2t) 



2t-2 4 2t-1 



+ 



;t-4 4)' 



where the coefHcients of the anomalous dimensions defined in Eq. (3.22) and the (3 function 
are in the MS scheme; they are given in Sec. 2.1 in Ref. [31]. 

Next, we apply the inverse MeUin transform of Eq. (3.19) in order to convert the 
resummed result to momentum space at NNLL accuracy (see Rcf. [31] and Sec. 3.4 in 
Ref. [36]). Let us consider the integrated width with r < tq for some tq < (1 — x)/\. The 
NNLL resummed result, matched to the full NLO expression takes the form: 



1 dT{X, r < ro, x) 



dXdx 

exp <! go{uj) 



[V{X,x) + AR'^{X,x)] X 

-1 



(3.26) 



1 + 



^g'MPi^^^^i^' (1 - Pog'oH) - (1 - PogM) ) 



TT 

'o 



TT 



where 



TT ro 



and where AR°°{X, x) is given explicitly in Eqs. (3.35) and (3.37) below. The regular term 

at 0{as) is included here directly in momentum space. 

Fig. 1 shows the differential spectrum as a function of r, for the chosen values of A = 0.9 
and X = 0.7, computed as a derivative of the resummed spectrum in Eq. (3.26). The large 
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Figure 1: The differential spectrum based on the fixed-logarithmic-accuracy formula of Eq. (3.26), 
which is matched to NLO. The LL, NLL and NNLL accuracy results are plotted as dotdashed, 
dashed and full line, respectively. The three curves end at r ~ 0.028, where the resummed results 
become complex owing to the Landau singularity at a; = 1/2. 



differences between curves corresponding to increasing logarithmic accuracy demonstrate 
the problematic nature of this expansion, in which infrared renormalons (in particular the 
one at u = ^) are unregulated. This problem is solved in what follows by computing the 
moment-space Sudakov exponent using Borel summation. 

Calculation of the Sudakov factor by Borel summation 

Direct evaluation of the Borel integral in the Sudakov exponent of Eq. (3.20) using the 
Principal- Value prescription makes optimal use of the known all-order structure of the 
exponent. In the context of DGE for inclusive spectra, this regularization has been proven 
effective in several respects [31]: 

• No Landau singularities are present. 

• It provides a systematic definition of the perturbative sum. Using this definition, 
cancellation of ambiguities can be realized — this has been explicitly used for the 
leading renormalon ambiguity associated with the pole mass. 

• The Sudakov factor is a real-valued function of the moment variable. 

• The resummed spectra has modified support properties. With the appropriate Bs{u) 
(see below), the support is close to that of the physical non-perturbative distribution. 
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The calculation of the Borel sum in Eq. (3.20) in a given regularization for the renor- 
malons, requires, in principle, the knowledge of Bs{u) and Bj{u) for any positive u. In 
contrast to the large-/3o limit (3.16), analytic expressions for these functions in QCD are 
not known. The fixed-logarithmic-accuracy approach described above uses only the expan- 
sion of these functions, as well as factors r(— 2-^) and T{—u) that multiply them, around 
the origin (tt = 0). As discussed in detail in Sec. 2.3 in Ref. [31], one may use additional 
information on the Borel functions Bs{u) and Bj{u) to constrain further the Sudakov fac- 
tor. Most importantly, Bs{u) at w = | was determined in Ref. [31] with a good accuracy 
using the cancellation [30] of the leading infrared renormalon ambiguity in the Sudakov 
factor with that of the pole mass. 

To understand the origin of renormalons in the Sudakov factor note that the integration 
over r near the r — > limit necessarily involves some contribution from the infrared region. 
^Prom Eq. (3.21) it is obvious that for small r the coupling is evaluated at small momentum 
scales, where it is out of perturbative control. In the Borel formulation, when computing 
the moments in Eq. (3.18) where the integrand is given by Eq. (3.17), the integration 
gives rise to the factors r(— 2^) and V{—u) in the soft and jet parts of the exponent, 
respectively. These factors contain infrared renormalons poles at integer and half integer 
values of u. This way infrared sensitivity translates into power-like ambiguity of the 
Borel sum [30,31,36 39], corresponding to powers of nA/{Xmh) and nA^/(Am5)^ in the 
soft and jet factors, respectively; the latter are obviously less important than the former. 
As shown in Refs. [30,31], the leading power ambiguity on the soft scale proportional to 
(n — l)A/(Am5), cancels exactly with kinematic power corrections. In the present context 
such corrections arise when expressing r in terms of the hadronic variables; see Sec. 4. 

Sublcading power ambiguities on the soft scale, corresponding to powers of [nK/ {\m},)y 
with^ j > 3, are related to the momentum distribution of the b quark in the meson. In 
principle, these ambiguities are removed only upon including a non-perturbative function, 
namely an infinite set of power corrections, that makes for the difference between the quark 
distribution in an on-shell heavy quark and that in a meson. The closer to the endpoint 
one cuts the distribution, the larger is the weight of high moments, and with it the sig- 
nificance of these power corrections. Our analysis of the B — > Xg'j decay [31] has shown 
that if Bs{u) is sufficiently constrained (see below) and kinematic power corrections are 
included, one can determine the photon-energy spectrum to a good accuracy neglecting 
any additional non-perturbative power corrections. 

In principle, the theoretical uncertainty associated with the unknown form of the 
anomalous dimension Bs{u) away from the origin, and the unknown non-perturbative 
power terms that characterize the quark distribution in the meson are two completely 
distinct issues. In reality, however, both result in similar, parametrically-enhanced power 
effects, and are therefore impossible to distinguish phenomenologically. 

The approach we take in this paper follows closely what was done in Ref. [31]: starting 
with the large-/3o result for Bs{u) and Bj{u), we include 0{u) and 0{u'^) terms to obtain 
the exact exponent to NNLO. For Bs{u) we further include O(u^) and higher-order terms 

'^The absence of ambiguity of 0{{nA/ (Xmb))^) is discussed below. 
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so as to match the computed residue of at « = i , getting^ 



Bs{u) = e3"(l - u) exp < C2U + 



1 



2^Ca /5 9 



^2 ^ X W{u), 



^ . , 1 5,, / 1 1 \ sinvru 
Bj(u) = -C3" \ X 

' 2 \1-U I-U/2J TTU 

43 



(3.27) 



2 / 29 2 — r/- \ , 



u 



2! 



+ 0(u^) , 



exp 102"^ + 

where 02,3 are given in Eq. (2.22) in Ref. [31] and where 

W{u) = e*i"+l*2«' ^1 _ +\{ti- ^2^^ = 1 + 0{u^). (3.28) 

Here ti,2 are fixed requiring: 

55(n= 1/2) = 0.91423 ± 3% (computed) ; ^5(^ = 3/2) = -0.23366 x C, (3.29) 

where we used the computed normaUzation of the leading renormalon residue of the pole 
mass at u = ^ (see eq. (2.30), (2.36) and (4.8) in [31] with Nj = 4) and an arbitrary 
normalization constant C at n = | . 

We stress that the ansatz of Eq. (3.27) assumes that Bs{u) vanishes at u = 1 and has 
no other zeros at positive u, as in the large-/3o limit (3.16). The constraint Bs{u = 1) = 
was found to be quite important [31]. Here we explicitly assume that it holds in the 
full theory and do not investigate the possibility of it being violated. The vanishing of 
Bs{u = 1) implies in particular that the perturbative Sudakov factor of Eq. (3.20) is free 
of 0{{nK/ {Xm,},))'^) ambiguities^. Thus, after the cancellation of the u = ^ ambiguity, the 
leading ambiguity is 0((nA/(Am;,))^), as already mentioned. To the extent that renor- 
malons do indeed give good indication of which non-pcrturbative effects are important, 
the absence of the u = 1 renormalon is well supported by the successful comparison of the 
prediction for the first two cut moments in B — > Xsj with experimental data. 

Note also that Eq. (3.27) tends to zero at asymptotically large u. In this respect it 
differs from the large-/3o limit. It is important to stress though that convergence of the 
Borel integrals at large u is guaranteed independently of this assumption, owing to the 
factors r(-2u) and r(-u) in Eq. (3.20). 

In Eq. (3.27) we parametrize Bs{u) at intermediate u values by a single parameter C. 
When Bs{u) is inserted into Eq. (3.20) this number determines the renormalon residue at 
n = |. The effect of C on the function Bs{u) is shown if Fig. 2. While the actual residue in 
the full theory is not known, it is expected to be smaller than in the large-/3o limit, where 
Bs{u = 3/2) ~ —6.09. Such a large residue can be obtained in the ansatz of Eq. (3.27) 



*This specific ansatz was called 'model (c)' in Ref. [31]. 

®A well-known analogous situation is the absence of the leading renormalon ambiguity [62] in Drell- 
Yan production [63]. Also there parametrically-enhanced power corrections appear, which are related to 
renormalon ambiguities in the Sudakov exponent [37]. However, based on the analytic structure of the soft 
Sudakov anomlaous dimension computed in the large-/?o limit, the first power correction is assumed not to 
appear. 
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Figure 2: The function Bs{u) with vertical Unes indicating the positions of possible renormalons 
owing to the poles of T{—2u) in Eq. (3.20). The thin line shows the large-/3o limit of Eq. (3.16) 
while the thick lines show Bs{u) according to Eq. (3.27) with Nf = 4 and with different assignments 
of the parameter C that controls the residue of the Borel singularity in the Sudakov exponent at 
M = |: C = 0.01 is shown as a dotdashed line while the default value C = 1 and — 2,4, 6, 8, 10, 12 
are shown as full lines. 



with iVj = 4 by setting C ~ 26 {W{3/2) ~ 20). In model (c) in Sec. 2.3 of [31], we had 
Bs{u = 3/2) ~ -0.23366 corresponding to C = 1 (TV(3/2) ~ 0.776). Here we keep C 
as a free parameter by which we gauge the theoretical uncertainty in the Sudakov factor. 
The possible range of values for C is discussed further in Sec. 4.2, where we also study the 
impact of this parameter on the spectrum and on the measurable partial width, see Fig. 8 
and Fig. 9, respectively. 



3.3 Matching and predictions for the logarithmic terms at NNLO 

Using the definition of the moments in Eq. (3.18) with the NLO expansion of Eq. (3.7) we 
have 



dTn{X,x) 
dXdx 



1 r"" 
— dr{l 

1 Jo 



dT{X,r,x) 



dr (1 — r) 



n-l 



dXdrdx 
V{X,x)6{r) + R{X,r,x) 



wo{X,x) + 



TT 



wi{X,x)+K^iX,x) 



where = (1 — x)/X denotes the maximal value of r and 



K^{X,x) 



dr {I — r)"' ^ki{X,r,x) 



(3.30) 



(3.31) 
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where we split the integration into singular and regular terms for r — > 0: 

Jo 



and 



(3.32) 



(3.33) 



rrm r 

= / dr (1 — r)"~ /ci(A,r, x) — fc^'"'^ (A, r, x) 
JO L 



where A;^'"*^ (A, r, x) and A;i(A,r, x) is given by Eq. (3.10) and Eq. (3.9), respectively. In 
Appendix B we evaluate if"^'°^ (A, x) and K""^ (A, x) explicitly. The result for the singular 
part reads 



K^^'^'i^x) = -wo{X,x] 



y-^'i(ra) + (5'(n) + 7£;)2 



-^(*(n) + 7E) (3.34) 



E 

fe=0 



(1 - r^)"+^ 
n + k 



ln(rm) + ^'(n + fe) - *(n) 



4 n + A; 

fe=o 



while the regular part is given in Eq. (B.13). 

In order to determine the n-independent coefficient that multiplies the Sudakov factor 
in Eq. (3.20) at ©(a^), 



Ai?-(A,x) = AK^{X,x)^^^^^^ + 0{al), 



TT 



(3.35) 



we compare the expansion of Eq. (3.20) with the NLO result of Eq. (3.30). We find that 
the 0{as) terms that are not generated by the expansion of the exponent are: 



AK^{X,x) = K^{X,x) 



Wo, 



(3.36) 



where L = In n, so 



AK^{X,x)= lim AK'^{X,x) = -woiX,x) 



12 ^ 2 4^^ 



(3.37) 



Using this coefficient Eq. (3.20) yields the following perturbative expansion: 



drn{X,x) 



dXdx 



wq{X,x) + 



large n 



7 



tI tt^ , 7^ 



12 4 



wq{X,x) 



-wi{X,x] 



CFas{mb) 



TT 



+ 



24 12 ' ' 



- 21 - 



+ 
+ 

+ 



IICa Nf 



12 6 

7E _ 13 \ 
4 144 

12 3 



Cf lnA + 

77 , 11 7^; 



21 3 2 49 7r2\ „ 

Y^^ + r^+32 + 24]^^ 



(3.38) 



24 6 



Cf InA 



,3 3 ^ 
+ 1^ + 2^3 



48 



+ 



16 



+ 



12 



+ 
+ 
+ 



85 5 2 13 1 2 \ 

\ TT '1E + -li) Nf 

144 72 72 ' 4 '^J ^ 
905 ^ ^ 957E _ 677r^ _ II7I _ 1 A c^)cp] 



288 12 
^ + 



144 



144 



7E + - ) CfL 



Cfwi{X, x] 



L 

TT J 



wo{X,x) 



All the coefficients of the logs determined here are exact. It is important to emphasize that 
virtual corrections at one-loop order, wi{X, x), contribute to log terms (CpL and C'^L?) at 
0{a^g). Similarly, two-loop virtual corrections, it;2(A,x), which are yet unknown, contribute 
to subleading logarithms at 0{a^g) and beyond. 

Complete matching at 0{as) can be achieved by including 0(l/n) terms additively, 
using "R matching" : 
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or, alternatively in front of the Sudakov factor, using "log-R matching" 



dTn(X, x) 



dXdx 



matched 



ViX,x) +AK-iX,x)^^^^+0{al] 



TT 



(3.40) 



,Cf f°°du^, , 
exp -f / —T{u) 



m^A^ 



[Bs{u)T{-2u) (n^" - l) - Bj{u)T{-u) (n" - 1)] 



/3o Jo u 

where the terms missing at 0{a^) include n-independent as well as 0{l/n) terms. 

Converting the expansion in Eq. (3.38) to momentum space we obtain for the log- 
enhanced contribution to the triple differential width (see Eq. (3.7)) at this order: 
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This expression is useful for checking NNLO calculations. An immediate application is the 
calculation of the log-enhanced part of partially integrated and single differential distribu- 
tions. For example, the single differential distribution with respect to p: 
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(3.42) 



The C(a^) coefficient that is leading in A^^ was recently computed [34] and it is consistent 
with Eq. (3.42). The other coefficients are new. 



3.4 Exponentiation beyond logarithms 

Let us note that the moment -space function obtained by integrating the r — > singular 
terms of the form r~^~^'^ or r~^~" in Eq. (3.17), which are regularized as ()* distributions, 
is not purely logarithmic. Performing the calculation leading to Eq. (3.20) above but 
avoiding any further larger approximation one obtains the following Sudakov factor: 



Sud(m&A,n) = exp 



(Of du 
\ Po Jo u 



T{u) 



A2 



m; 



'A2 



(3.43) 



Bs{u)T{-2u) 



Tin) 



1 



r(n - 2u) r(l - 2u) 



-Bj{u)ri-u) 



r(n) 

r(n — I 



r(i 



-23- 



Note that to any order in u: T(n)/T{n — u) ~ n" up to constant terms. Thus, upon 
applying this approximation one returns to Eq. (3.20). There, only powers of Inn appear 
at any given order in u, while constants and 0{l/n) terms are entirely excluded from the 
exponent and appear only in the matching coefficient. In Eq. (3.43) such terms appear also 
in the exponent. 

It is straightforward to match the Sudakov factor of Eq. (3.43) to NLO. The result 
is more elegant than when using a "purely logarithmic" Sudakov factor (Eq. (3.40) or 
Eq. (3.39)) because the constants for n — > oo, which where included there in the process 
of matching, appear in Eq. (3.43) as part the exponent. This is true at higher orders too: 
when using Eq. (3.43), AR°° = 0. Consequently, it is just the purely virtual terms ^(A, x) 
(terms proportional to 5{r), after the ()+ distributions have been defined) that must be 
multiplied by the Sudakov factor in order to generate correctly the log-enhanced terms at 
higher orders. Using Eq. (3.43) with "R matching" we get 
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while with "log-R matching" we get 
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where the last line corresponds to K^^' obtained in Eq. (B.13). The definitions of the 
coefficients fi{X,x) and fi{X,x) and the integrals In+i{rm) and Jn+i{rm) are summarized 
in Appendix B. 

Recall that all the terms in AK^{X,x) are suppressed by at least one power of 1/n. 
Note, on the other hand, that for fixed n this matching term contains logarithms of r^, 
which get large for — > 0. In particular, the first moment (n = 1) is: 



n=l 

AK^ {X,x) = wq{X,x) 



1 7 
-\v?{rm) - -ln(r„,) 



+ Kr^"«(A,x), (3.47) 



where i^J""' is given in Eq. (B.9). 

It is interesting to study the renormalon structure of the Sudakov factor in Eq. (3.43) 
and compare it to the simplified version of Eq. (3.20). Let us examine in some detail the 
soft function (a similar structure appears on the jet side). As discussed in Sec. 4.1 in 
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Ref. [31], power corrections associated with renormalon ambiguities in the soft function 
modify the Sudakov factor of Eq. (3.20) or Eq. (3.43) as follows^'': 



Sud(m6A, n)|j 



Sud(m6A,n)|py x 



(3.48) 
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where the n-dependent residue function in the two cases is 
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and fj^ are dimensionless non-perturbative parameters of order 1. PV stands for the 
Principal Vahies prescription, which is used to define the perturbative sum on the one hand 
and the non-perturbative parameters on the other. In Eq. (3.48) the l.h.s. is regularization 
prescription dependent, while the r.h.s. is not. 

It is important to note that subleading renormalons O (j^A/{Xmb)y^ corresponding 
to increasing j have an additional numerical suppression owing to the structure of the 
Sudakov exponent, namely the fact that 



Resr(-2«)|„^./2 = - 



1 (-l)-'+i 
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Considering only powers of nA/{\nn,) the residue structure of the two exponents is the 
same. In order to analyze the differences between them at smaller n let us now compare 
between the residue function R{n,j/2) in the two cases for the first few renormalons at 
u = j/2: 
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(3.51) 



Obviously, using Eq. (3.43) the parametric enhancement of power corrections is not as 
dramatic as it looks based on Eq. (3.20). In Eq. (3.43) integer moments are entirely free 
of certain power-like ambiguities which do show up in Eq. (3.20): according to the residue 
structure of Eq. (3.43), for n = N, N being a positive integer, only powers {A/{Xmb)y 
with j < N appear. 



'^"As discussed in Sec. 4 below, the j = 1 ambiguity cancels upon converting to hadronic variables and the 
j = 2 renormalon is absent since -65(1) = 0. Thus, dynamical non-perturbative power corrections actually 
appears only for j >3. 
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A priori, one might expect that ah power corrections would become of order one at 
n ~ Xrrib/A. Assuming that renormalon ambiguities do indeed give good indication on 
non-perturbative effects, the factorial suppression of Eq. (3.50) as well as the n depen- 
dence of the residues of Eq. (3.43) imply that power terms are altogether less important 
and, in contrast with the naive expectation, the power expansion corresponding to higher 
renormalons in Eq. (3.48) does not break down. 

3.5 Resummed double— differential width with a lepton— energy cut 

Since the dependence on the lepton energy fraction x appears in the resummation formula 
(3.45) only though the phase-space limits and the matching coefficients but not through 
the Sudakov factor, it is possible to integrate over x analytically to obtain the double 
differential width with a cut x > xq: 
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In the Appendix C we explicitly compute the matching coefficients at ©(a^), namely, 
Vr{X, X > Xo) and AKi^{X, x > xo), by integrating the virtual coefficients of Eq. (3.7) and 
the real-emission moments of Eq. (3.46), respectively. 



4. Pctrtial width with experimentally relevant cuts in hadronic Vciriables 
4.1 Conversion to hadronic variables and power corrections 

So far we examined the resummed distribution in kinematic variables associated with the 
partonic jet initiated by the u-quark. Measurements do not distinguish between the u-quark 
jet and the soft partons originating in the light degrees of freedom in the meson. Working 
still in the approximation where the b quark is on shell, one defines v = pb/Mb (where 
= 1) so = fTT-hV and the momentum of the light degrees of freedom is At; = (M^ — m^)!;. 
The hadronic variables (in the B rest frame) are then: Eh = Ej + A, where A = Mb — mh, 
and Ph = Pj- The relations with the corresponding hadronic lightcone variables are {cf. 
Eq. (3.2)): 



P+ = Eh ~\Ph\ = A + Ej - \pj\ = mbp + A, 
P~ = Eh + \Ph\ = A + Ej + \pj\ = rubX + A. 



(4.1) 
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In the small hadronic-mass-squared {P~^P~) region, p is small, and the A term in P"*" is 
absolutely essential. On the other hand, the subregion where also A is small {P^ is of order 
A) is unimportant because it is power suppressed by the Born-level weight (see Eq. (C.2)). 
Nevertheless, as we carefully treated all other 0{A/mi,) effects, we do not neglect A com- 
pared to m^A when converting the decay width to hadronic variables. 
In terms of hadronic variables the total phase space is: 



Jo Jmb-2Ei Jo dP+dP dEi 



(4.2) 



A perturbative calculation of the differential width obviously cannot fill the entire hadronic 
phase space. The fixed-order perturbative result simply has a smaller phase space, c/. 
Eq. (3.4). However, resummation can, and indeed does modify the support properties, 
violating energy and momentum conservation which hold order by order. In this respect 
it is essential to distinguish between the region targeted by the resummation and other 
phase-space regions: 

• In the small P+ region (or, equivalently, the small hadronic mass region) the resum- 
mation-*--^ dramatically improves the description of the differential width, and with 
it the support properties. The resummed perturbative distribution of Eq. (3.19) 
extends to negative r values. In terms of hadronic variables, the differential width 
extends to the region P+ < A and approximately vanishes for P+ < (see Fig. 3), 
consistently with the physical support properties [31]. 

• Away from the region targeted by the resummation, partonic phase-space boundaries 
are also violated^^. This, however, is an artifact of uncontrolled 0{l/n) higher-order 
perturbative corrections, which depend on the specific scheme by which matching 
to the fixed-order result is done. As shown in Fig. 3 the phase-space limit r < 
(1 — x)/X (see Eq. (3.4)) is violated. Fortunately, this makes just a small impact on 
the observables we consider. This issue is discussed further in Sec. 4.2. 

Considering other phase-space limits of a; = 2Ei/mi, or A = (P~ — A)/mb that are 
not related to the resummation variable r, the partonic phase-space limits of Eq. (3.4) are 
unmodified by the resummation. The phase space after the P"*" integration is shown in 
Fig. 4, where we indicate both the hadronic and the partonic upper limits^^ on the charged 
lepton energy. 

As explained below, the transformation of the resummed distribution in the small 
lightcone component from partonic to hadronic variables has a crucial role in obtaining a 



'^^Note that we refer here specifically to DGE where renormalons are regulated by the Principal Value 
prescription. This improvement is not achieved in fixed-logarithmic-accuracy Sudakov resummation, see 
e.g. Fig. 1, which suffers from Landau singularities. 

Violation of the hard phase-space boundary by resummed distributions is a familiar phenomenon, see 
e.g. Rcf. [38] 

^■^Note that although we shall refer below to the hadronic phase space, Ei < Mb/2, the numerical 
integration is actually done in the range Ei < mb/2, as dictated by the perturbative support properties. 
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Figure 3: The fully differential width (in units of Tq) computed based on Eq. (3.45), plotted as 
a function of the small lightcone component P+ at P~ = 3.5 GeV and at three different values of 
El = 1, 1.5 and 2 GcV. The maximal P+ endpoint {Mb — "^Ei), denoted by a vertical line, should 
be: 3.3, 2.3 and 1.3 GeV, respectively. The u = 3/2 renormalon residue parameter in Eq. (3.27) (see 
Fig. 2) is C = 1. Note that the perturbative support properties, namely A < P"*" < Mb — ^Ei are 
violated at both ends. At the low P"*" end this is of course a desired consequence of the resummation: 
the physical support properties P+ > are approximately recovered. At the high P+ end this is 
an 0{al/n) artifact, which is non-negligible only for large Ei. 



correct, unambiguous result (which is also characterized by approximate physical support 
properties). 

To this end, let us recall some of the observations of Refs. [30,39] and reformulate them 
in the context of the semileptonic decay. The first observation, already discussed above, 
is that renormalons in the soft function of Eq. (3.20) have residues proportional to inverse 
integer powers of the soft scale, indicating potential non-perturbative power corrections 
of the form (nA/(Xmh)y, where j is a positive integer. It was also observed there that 
the dominant part of these corrections should be attributed to kinematic effects associated 
with the mass difference between the meson and the quark, A = Mb — mi,. In the present 
context these correspond to identifying r = Pj' /pj = p/A in terms of hadronic variables 
using Eq. (4.1): 

P+ -A , , 

which effectively shifts the distribution in r by A/{Xmf,). Note that Eq. (4.1) relies on the 
assumption that the b quark is on-shell, so this transformation is purely kinematic and 
the "primordial" motion of the b quark in the meson is not taken into account. In order 
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Figure 4: The full phase space in the Ei) plain (upper left triangle) after the P+ integration 
has been performed. The dotted horizontal line Ei — mh/2 starting at Pj^ = A is the upper 
boundary of the partonic phase space. A lower cut on the charged lepton energy (here Ei > 
i?o = 1 GeV) removes the small triangular region below the horizontal dashed line. An upper 
cut on the hadronic invariant mass (here P+P^ < = (1.7 GeV)^) affects the distribution for 
El < ^ {Mb — M"^ I P~ ) . A full line separates between this region and the rest of the phase space 
where this cut is irrelevant. On this line we denoted three special points by D [P~ = Mx), C 
{p- = Mb - 2Eo) and F {p- = Mb), respectively. 




to account for such non-perturbative dynamics, additional power corrections of the form 
{nA/{Xmb)y with j >2 should be included, as discussed below. 

Still within the on-shell heavy-quark approximation, converting Eq. (3.19) to hadronic 
variables using Eq. (4.1) one obtains: 



1 dT{P+, P-, El) 
To dP+dP-dEi 
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dXdx 



ml P 



C+lOO 

c—ioo ^"^ 
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dTn{X 



mi, 



2E, 



dXdx 



(4.4) 



where no further approximation was made. The A dependence associated with the shift 
in P~^ factorizes exactly in these variables. Note, in particular, that the kinematic power 
correction factor (l — A/P")" is valid for any n, not just for large n. Obviously, it makes 
no impact of the first moment and it becomes increasingly important as n gets large. 
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It was further shown^^ in Ref. [30] that the product 



d\dx \ P~ ) d\dx 

is free of the leading renormalon ambiguity, a hnear ambiguity of O (A(n — 1)/P~) in the 
exponent, owing to exact cancellation of the u = ^ ambiguities between the Sudakov factor 
in the resummed moments and the one in A. The exponentiation of the corresponding 
power-like ambiguity in Eq. (3.20) and the regularization of A and of the Sudakov exponent 
using the same prescription are crucial for this cancellation to take place. In the Principal 
Value prescription the numerical result for A is 

A = Apv = Ms - ml"" = 5.28 - 4.88 ± 0.06 GeV = 0.40 ± 0.06 GeV, (4.6) 

where we used the result for the mass ratio in Eq. (2.4) with the short-distance quark mass 
value of Eq. (2.11) and Nf = 4. 

Let us address now power corrections associated with the dynamical structure of the 
meson. The perturbative calculation of the soft Sudakov exponent accounts for radiation off 
the heavy quark that puts it slightly off its mass shell. However, it cannot take into account 
the way in which the virtuality of the heavy quark is influenced by its non-perturbative 
interaction with the light degrees of freedom in the meson. As discussed in Refs. [30,31], in 
the present framework such non-perturbative dynamics is reflected in additional power cor- 
rections of form (nA/(Xmb)y with j > 2. The detailed structure of these power corrections 
can be read off from the renormalon ambiguities of the soft Sudakov exponent. As done 
explicitly in Eq. (3.48) each power ambiguity is associated with a new non-perturbative 
parameter fj^; upon including the power term, the corresponding ambiguity is removed. 
The multiplicative correction to the Sudakov factor in Eq. (3.48) translates directly into a 
multiplicative correction to partonic moments computed, for example, in Eq. (3.45). Given 
the kinematic cancellation of the leading (u =1/2) renormalon ambiguity discussed above, 
and given that B^iu = 1) vanishes (see Eq. (3.27) and the discussion following it), the lead- 
ing dynamical power correction is associated with the u = 3/2 ambiguity, corresponding 
to j = 3. This correction would modify the spectral moments by 



drn{x,x) 

dXdx 



(ir„(A,x) 



dXdx 



X (4.7) 



PV 




where R{n, 3/2) = ^(n — 3)(n — 2)(n — 1). Note that unless the residue of the anomalous 
dimension 5^(3/2) and the non-perturbative coefficient /J"^ are much larger than one, 
which we consider unlikely, Eq. (4.7) amounts to a small correction even for n ~ Am{,/A. 

The corrections of Eq. (4.7) go beyond our minimal model: here we regularize all 
the renormalons using the Principal Value prescription and do not include any additional 

^''The cancellation has been explicitly checked by a calculation in the large-/3o limit and argued to be 
general. 
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non-perturbative power terms (i.e. we set fj^ = for j > 3). The resulting, essentially 
perturbative, spectrum of Eq. (4.4) may therefore differ from the measured spectrum by 
effects of the from of Eq. (4.7); because of the parametrically small contribution to moments 
n ^ Xrub/A, these effects are restricted to a narrow region near P+ 2± 0. We shall revisit 
this issue in Sec. 4.2 below when estimating the related theoretical uncertainty. 

In Appendix D we explicitly show how the resummed result for the spectrum, computed 
in the previous section using moments of the partonic lightcone momentum ratio r = p/X, 
is used to compute the partial branching fraction with experimentally-relevant cuts. In 
Sec. D.l we use the results of Sections 3.3 and 3.4 to express the partially integrated 
width with a cut on Mx (or on P+) and an additional mild cut on the lepton energy. In 
Sec. D.2 we derive expressions for the same observables using the results of Sec. 3.5 and 
Appendix C, where the cut on the lepton energy is implemented analytically. Numerical 
results for i?cut as a function on Mx or arc presented in Sec. 4.2, where we also perform 
a detailed study of the theoretical uncertainty. Finally, in Sec. 4.3 we extract \Vub\ from 
recent measurements by Belle. 

4.2 Numerical results and theoretical uncertainty estimates 

In extracting \Vuh\ from experimental data according to Eq. (1.2) the effect of kinematic 
cuts is taken into account through the event fraction i?cut defined in Eq. (1.1). Here 
we shall use^^ the tools of the previous sections to evaluate Rcut for hadronic mass Mx or 
lightcone momentum P"*" cuts (with additional lower cut on Ei) and estimate the theoretical 
uncertainties involved. 

It should be emphasized that although the formulae in Sections D.l and D.2 are ex- 
plicitly written for partial rates and are normalized by Tq, we are really dealing here with 
the perturbative expansion of Pcut^ not the partial width and the full width separately. 
A single calculation of Pcut in Eq. (1.1) involves using these formulae twice: with a cut 
in the numerator, and with no cut in the denominator. Clearly, Fq cancels in the ratio. 
This has far-reaching implications in what concerns renormalization-scale dependence and 
renormalon ambiguities: Pcut, in contrast with the partial (or total) rate in units of Fq, is 
not affected by the ©(A/mt) renormalon ambiguity. This ambiguity would have cancelled 
with the ambiguity of in Fq had this factor been included and the series resummed. In 
Pcut it simply cancels^^ between the numerator and the denominator. 

As already mentioned, we focus on the cut that maximizes the rate, as measured by 
Belle [7]: P+P~ < M^ = (1.7 GeV)^ to discriminate the charm background with an 
additional lepton-energy cut Ei > Eq = 1 GcV. Numerical integration of the differen- 
tial distribution with respect to P~ and Eq. (D.2), over the available phase space 
yields Pcut(-^x = 1.7GcV,£'o = 1 GeV) = 0.615. This event fraction will be used in the 
next section to extract from Belle data. In this calculation we used the matching 

of Eq. (3.45) where the renormalization scale for the coupling was set as /it = P~ . In the 

Numerical analysis is done using a C program that combines Borel integration, Mellin inversion and 
phase-space integration [64]. 

^^Note, however, that cut-related renormalon ambiguities, which are parametrically enhanced at large n, 
are not canceled in this way; these are discussed below. 
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soft Sudakov exponent we chose C = 1. In the following we perform a detailed study 
of the theoretical uncertainty in Rent- We consider uncertainties from three sources: (1) 
higher-order perturbative corrections going beyond the resummation as well as beyond the 
NLO; (2) renormalons and power corrections, where the main uncertainty is associate with 
parametrically-enhanced power terms related to the quark distribution in the meson; (3) 
uncertainty in the values of fundamental short-distance parameters: af^ and m^^. 

Matching schemes, scale dependence and higher— order perturbative effects 

The resummation employed here focuses on improving the perturbative expansion of the 
fully differential width in the particular kinematic region where one lightcone component of 
the hadronic system, P+, is small. Matching has been performed to the full NLO expres- 
sion. Our purpose here is to estimate the effect of higher-order perturbative corrections 
on i2cut- This is done in two ways: first by comparing different matching schemes in- 
troduced in sections 3.3 through 3.5, and then by renormalization-scale variation in the 
matching coefficient. 

Let us first estimate the numerical significance of the difference between different 
matching schemes, which are formally of 0{a^/n). Fig. 5 shows the result of three differ- 
ent calcTilations of the event fraction .Rcut with a hadronic mass cut P~^P~ < M|- as a 
function of Mx, with a fixed lepton-energy cut Ei > 1 GeV, all with the same choice of 
renormalization scale for the coupling in the matching coefficient, /j, = P~ (see below): 

(a) Using the partially integrated distribution with respect to Ei, namely Eq. (D.ll) and 

Eq. (3.53). 

(b) Using the differential distribution with respect to Ei, Eq. (D.2) (or, equivalently, Eqs. 

(D.4) through (D.6)) with the matching of Eq. (3.45). 

(c) same as (b) with the matching of Eq. (3.40). 

The origin of the difference between (c) and (b) is clear: the constant terms, which were 
treated in Eq. (3.40) as part of the matching coefficient are exponentiated in Eq. (3.45). Less 
obvious is the difference between (a) and (b), which is reflected also in the P" distributions 
plotted in Fig. 6. Let us explain this difference for the example where an Mx cut is applied. 

The first crucial observation is that the resummed distribution (Fig. 3) does not respect 
the phase space limit r < {1 — x)/X, which holds order by order in perturbation theory. 
Instead, the distribution develops a tail that extends up to r = 1. By construction, the first 
moment, which is determined by the matching coefficient (i.e. the fixed-order result) alone, 
is unaffected by this tail. In contrast, whenever one uses the distribution computed as an 
inverse-Mellin of the resummed moment-space expression these details become important. 

Considering next the way the phase space in Fig. 4 is covered we see that in case (a) 
only the region P~ < Mx is evaluated using the first moment whereas the entire region 
P~ > Mx is reconstructed as an inverse Mellin transform of the resummed expression. 
In case (b) the differential distribution is used only where necessary: the entire region 
above the line Ei > \{Mb — M\/P^) is computed using the first moment and only the 



-32- 




Figure 5: The event fraction i?cut within the range P^P^ < AI^ and Ei > 1 GeV, plotted 
as a function of the cut value Mx- Three different matching schemes are used: (a) the partially 
integrated distribution with respect to Ei, namely Eq. (D.ll) and Eq. (3.53), numerically integrated 
over (full line); (b) the fully differential distribution computed by Eq. (D.2) and the matching 
of Eq. (3.45), numerically integrated over the Ei and P~ (dashed line); and (c) same as in (b) but 
with the matching of Eq. (3.40) (dotdashed line). 



matching scheme 


Event Fraction Rent {Mx = 1.7 GeV, Eq = I GeV) 


variation(%) 


(a) 


0.5941 


-3.5 


(b) 


0.6154 


default 


(c) 


0.6364 


+3.4 



Table 1: Results of the event fraction i?cut within the range Mx < 1.7 GcV and Ei > 1 GeV using 
different matching schemes (see text) with the scale of the coupling in the matching coefhcient set 
as /i = P^ . 



region below this line depends on the details of the differential distribution. Consequently, 
(a) effectively uses the distribution (computed as an inverse Mellin transformation) in the 
region 

Mx <P- < min |mb, ^| , Ei > ^ (^Mb - ^ 

where (b) uses the n = 1 moment instead. 

Table 1 summarizes the results corresponding to P~^P~ < = (1.7 GeV)^. We see 
that for the P~^P~ < = (1.7 GeV)^ cut, the matching-scheme dependence amounts to 
less than ±4% uncertainty. Although the fully differential calculation (b) is used to obtain 
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Figure 6: The distribution after integration over P+ (using the moment-space result) and 
El in four different situations: no cuts (full line), Ei > 1 GeV cut only (dotdashed line), Ei > 1 
GeV combined with P+P^ < Mj. = (1.7 GeV)^ cut (dashed), and finally Ei > 1 GeV combined 
with P+ < P,Jax — 0-66 GeV cut (dotted). In each of the two last cases, two different calculation 
procedures are used in order to gauge the sensitivity to higher-order corrections (matching-scheme 
dependence): (a) using the partially integrated distribution with respect to Ei, namely Eq. (D.ll) 
and Eq. (3.53) and (b) numerically integrating over the differential distribution with respect to Ei, 
Eq. (D.2) and Eq. (3.45). As explained in the text, the difference between (a) and (b) exists only 
when there is a Mx or a Pj^^x cut and only for P~ > Mx, or P~ > PjJaxJ respectively. 



the central value, we will estimate other sources of uncertainty using scheme (a), in which 
numerical results are obtained faster. 

Next, let us consider the renormalization scale dependence. When considering the par- 
tial width with stringent cuts, Sudakov resummation is essential. As shown in Refs. [30,31] 
and Sec. 3.2 above, the logarithmic accuracy criterion fails because of large running-coupling 
effects. By using DGE we resum these effects as well. However, this resummation is re- 
stricted to the Sudakov exponent. A priori, one might worry that additional resummation 
of running-coupling effects should be applied in the matching coefficient as well. Using 
scale variation we show below that this resummation is probably not necessary for the 
calculation of i?cut at the ±5% level. 

When considering i?cut with very mild cuts Sudakov resummation is irrelevant. But 
running-coupling effects, which we have not resummed, may still be significant. It is 
important to emphasize, that these effects cannot be estimated using the expansion of the 
total width^^, Eq. (2.1), which in contrast with -Rcut; is affected by the infrared renormalon 

^'^Indeed, the BLM scale of the total width is of order A, see e.g. [56]. 
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Figure 7: The event fraction i?cut within the range P+P^ < and Ei > 1 GeV, computed 
using Eq. (D.ll) and Eq. (3.53), plotted as a function of the cut value Mx- Different assignments 
for the argument of the coupling in the matching coefhcient are used. 

corresponding to m^. Determination of the BLM scale in Rent requires the knowledge of 
the matching coefficients of Sec. 3.3 as a function of A and x to O(a^/?o)- These calculations 
have not yet been performed. 

We therefore resort to estimating running-coupling effects in the matching coefficients 
by means of scale variation. A natural scale for the coupling in the matching coefficient is 
P~ , which is the hard scale characterizing the hadronic system, /i = will be our default 
value. Clearly, as good an alternative would be the partonic pj — this scale appear as the 
hard scale in Sudakov exponent, see e.g. Eq. (3.20) or Eq. (3.24). 

Fig. 7 presents the result for the event fraction Rent associated with a hadronic mass 
cut P+P- < M\ as a function of Mx with a fixed lepton-energy cut Ei > 1 GeV, 
computed according to Eq. (D.ll) and Eq. (3.53), with three different assignments of the 
renormalization scale ^ in a^'^(/i) in both the virtual and real-emission contributions to 
the matching coefficient. The results for P+P~ < Mj^ = (1.7 GeV)^ (and Ei > I GeV) are 
also summarized in Table 2. We find that while the coupling varies from as(m{,) ~ 0.22 to 
as{mb/2) ~ 0.29, i.e. by over 32%, i?cut varies only by about 8.3%. 

Let us consider now the default choice of ^ = P" . As shown in Fig. 6 the P^ 
distribution in which no Mx (or P+) cuts are applied is peaked close to P^ ~ Mb- On 
the other hand the Mx cut suppresses the distribution at the hard-P~ end, while it makes 
no impact on the soft end. Therefore, the choice of matching-coefficient scale as /i = P~ 
is a priori quite different from = rrih- Nevertheless, as shown in Table 2 the overall effect 
of this change on i?cut for P+P- < M^ = (1.7 GeV)^ (and Ei > 1 GeV) is just ~ 3.9%. 
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scale of as 


Event Fraction R^nt {Mx = 1-7 GeV, Eq = I GeV) 


variation (%) 


fi = p- 
tJ, = mb 
H = mb/2 


0.5941 
0.6173 
0.5677 


default 
+3.9 
-4.4 



Table 2: Results of the event fraction within the range P+ < = (1.7 GeV) ^ and Ei > Eq ^ 
1 GeV using Eq. (D.ll) and Eq. (3.53) with different assignments for the scale of the coupling in 
the matching coefficient (the n = 1 moment). 

Taking into account the scale uncertainty as ±4.4% and combining it in quadrature^* 
with the matching scheme dependence estimate of ±3.5%, we assign a total theoretical 
uncertainty of 5.6% in i2cut for Mx = 1.7 GeV owing to unknown higher-order corrections. 
We expect these number would significantly reduce once NNLO calculations of the fully 
differential width would be completed. 

Renormalons and powder corrections 

Considering the leading, u = 1/2 renormalon ambiguity, it is essential to distinguish be- 
tween 

• the 0{A/mb) renormalon in the total width when written in units of Fq (or the n = 1 
moment of the matching coefficient) which cancels by defining i?cut as we did, and 

• the cut-related leading 0{{n — l)A/mb) renormalon in the Sudakov exponent which 
makes all the perturbative moments (but n = 1) ambiguous and generates an am- 
biguous shift of the normalized distribution. 

While the former is canceled in Rent by definition, the latter affects only the partial width 
(i.e. the numerator in -Rcut) and becomes more significant the deeper the cut is. As 
discussed in detail in Sec. 4.1, its cancellation [30,31] requires renormalon resummation 
in the Sudakov exponent and incorporating the corresponding kinematic power corrections 
(l — A/P")" ^ according to Eq. (4.4). Having full control of the u = 1/2 renormalon, 
including the normalization of its residue, and having defined both the Sudakov factor and 
A using the Principal Value prescription, the computed spectrum is free of any 0{A/ (Am^)) 
artifact of the non-physical on-shell heavy-quark state. 

Beyond the u = 1/2 renormalon the Sudakov factor has subleading power corrections 
reflecting the non-perturbative dynamics of the b quark in the B meson. As discussed 
in Refs. [30,31] and in Sections 3.4 and 4.1 above, these corrections depend on the definition 
of the perturbative sum applied in the Sudakov exponent. Moreover, the perturbative 
Sudakov exponent (3.45) as well as the corresponding renormalon ambiguities (4.7) depends 
on the magnitude of the Borel function away from the origin, so the assumptions made on 
Bs{u) are directly relevant. Let us briefly summarize these assumptions: 

'^^By regarding these two uncertainties as independent we obtain a conservative uncertainty estimate. 
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(1) We assumed that, as in the large-/3o hmit (3.16), the Borel function Bs{u) vanishes 
at M = 1 and has no additional zeros. The vanishing of -85(1) implies, in particular, 
that there is no ambiguity 0{{nA/ (Amfe))^) in the exponent so the leading correction, 
corresponding to the u = 3/2 renormalon ambiguity, takes the form of Eq. (4.7); 

(2) We assumed that Bs{u) gets small at large u, although this is not required for the 
convergence of the Borel integral. To control the contribution from intermediate 
values of u, we included in the ansatz for Bs{u) in Eq. (3.27) a single free parameter 
C which is proportional to Bs{u = 3/2); see Fig. 2 and Eq. (3.29). It is difficult to 
compute Bs{u = 3/2), but numerically large (as well as extremely small) values of C 
can probably be excluded by the following considerations: (a) for C ~ 26 the residue 
at u = 3/2 would be as large as in the large-/3o, which is unlikely since non-Abelian 
corrections tend to reduce renormalon residues; (b) given the ansatz of Eq. (3.29), 
either large or small values of C imply that ti and t2 become unnaturally large, 
and so do higher-order perturbative corrections to Bs{u). Eventually, these can be 
computed explicitly to further constrain this function. 

Having fixed the Borel transform of the anomalous dimensions in Eq. (3.27), we compute 
the Sudakov exponent of Eq. (3.43) applying the Principal Value prescription to all the 
Borel singularities. The resulting DGE spectrum displays only small sensitivity to C, see 
discussion below. Therefore, despite the formal dependence on the regularization prescrip- 
tion for the u = 3/2 and higher renormalons, this resummed spectrum can be directly 
considered an approximation to the meson decay spectrum. On the theoretical level, this 
is supported by the following observations [31]: 

• According to the renormalon ambiguity pattern, the leading correction corresponds 
to the third power of nA/{\mi,). Eq. (4.7) implies that for not-too-high moments 
{n ^ Xirib/A) the non-perturbative effect is small, and therefore only a narrow region 
near P"*" ~ is affected. 

• Despite the finite gap between the physical support properties (P"*" > 0) and the 
perturbative ones (P^ > A), the resummed perturbative spectrum, quite remarkably, 
has support that is close to the physical spectrum. 

The final judgement, however, must be based on comparison with experimental data. Com- 
paring the predictions of Ref. [31], corresponding to C = 1 in Eq. (3.27) above, with exper- 
imental data for B — > Xs'J, shows [41] that this minimal model, which essentially has no 
free parameters, describes the spectrum well. In particular, the first two central moments 
with cuts, namely the average energy and the width defined with Ej > Eq (Eqs. (5.3) 
and (5.4) in Ref. [31]) agree with experiment within error, while for the third moment the 
data is still statistically limited [12,13]. It should be emphasized that the theoretical error 
estimate in [31,41] did not include any variation of Bs{u) nor any non-perturbative power 
corrections, but only variation of the short-distance parameters a^'^ and m^^ within their 
ranges of uncertainty. 
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Figure 8: The fully differential width (in units of Tq) computed based on Eq. (3.45), plotted as 
a function of the small lightcone component -P+ at Ei = 1.5 GeV and P~ = 3.5 GeV, for three 
different assignments of the u = 3/2 renormalon residue parameter C in Eq. (3.27) (see Fig. 2). 
C = 1 is the value used in Ref. [31]. 



Given the good description of the B — > Xg'j spectrum [41] in the present framework, 
our default choice here is not to include any parametrization of additional non-perturbative 
corrections of 0{{nA/{Xmb)y) with j > 3. As soon as theoretical predictions and experi- 
mental results become sufficiently constraining, parametrization of these corrections along 
the lines of Eq. (4.7) would be in place. 

We note that since the size of the leading correction in Eq. (4.7) is controlled by 
the product of f^^ and Bs{2>/2) (or, equivalently, C) these two parameters are strongly 
correlated. Based on data alone it would therefore not be possible to distinguish between 
genuine non-perturbative effects in the meson (/J^) and unknown higher-order corrections 
to Bs{u) that modify its behavior away from the origin with respect to Eq. (3.27). To 
determine /J^, one would need not only very accurate data but also further theoretical 
constraints. 

Let us therefore proceed by analyzing the numerical consequences of varying Bg (u) in 
Eq. (3.43) by means of C, while setting fj^ (j > 3) in Eq. (4.7) to zero. Later we shall 
consider the uncertainty associated with a non-zero /J^. C influences the width of the 
spectrum as a function of P^, but even more, the left-right asymmetry in the shape. This 
can be seen in Fig. 8, where it is also obvious that large values of C, such as C = 10, are 
disfavored as the distribution extends then further into the negative domain, violating 
the physical support properties. 

The final effect of varying C within the range 0.1 to 10 on Rcut is shown in Fig. 9 
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c 


Event Fraction i?cut {Mx = 1-7 GeV, = I GeV) 


variation (%) 


0.1 


0.5799 


-2.4 


1 


0.5941 


default 


10 


0.5870 


-1.2 



Table 3: Results of the event fraction within the range P+ < Mj. = (1.7 GeV)^ GeV and 
Ei>l GeV using Eq. (D.ll) and Eq. (3.53) with different assignments of C in Eq. (3.29). 
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Figure 9: The event fraction i?cut within the range P'^P~ < M|- and Ei > 1 GeV, computed 
using Eq. (D.ll) and Eq. (3.53), plotted as a function of the cut value Mx- Different assignments 
for the unknown parameter C, which controls the shape of the perturbative quark distribution 
function (see Eq. (3.29)) are used. In the matching coefRcient the renormalization scale ^ = P^ is 
chosen. 

as a function of Mx- Results for the experimentally-relevant value of P^P^ < M|- = 
(1.7 GeV)^ (and Ei > Eq = 1 GeV) are also summarized in Table 3. We find that varying 
C from 0.1 to 10 the event fraction varies by only 2.4%. Note that the variation is non- 
monotonous. 

Finally, let us return to the numerical effect of the leading non-pertrbative correction 
in Eq. (4.7). As discussed above, one expects fj^ to be of 0(1). For /J^ = 1 and for our 
default value C = 1 the effect is very small. For example, for the experimentally-relevant 
cut value of _P+P" < Mj. = (1.7 GeV)^ the correction to i?cut is ~ 0.2%; naturally, 
the effect increases significantly when lowering the cut, for example, for P^P~ < M'j^ = 
(1.4 GeV)^ it amounts to 1.3%. As obvious from Eq. (4.7), if one assumes C = 10 instead, 
the corrections are larger. In this case with /g ^ = 1 one obtains 0.8% and 4.0% for the two 
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parameter 


Hjvent rraction .rLcut y^^^x = i-' *^ev, £/o = i i^ev j 


variation (%) 


MS ( \ n 1 1 c 

al [raz] = O.llo 


O.oooi 


— 1.8 


Ois [mz) = 0.119 


0.5941 


deiault 


„,MS /™ ^\ n 1 00 

ctg ymz) — u.izz 


u.oyoD 


U.o 


iV/ = 3 


0.6008 


+1.1 


Nf = A 


0.5941 


default 


^MS ^ 4 


0.5511 


-7.2 


mp = 4.19 


0.5941 


default 


mp = 4.24 


0.6325 


+6.5 



Table 4: Results of the event fraction within the range P+P < = (1.7 GeV)^ and -E; > 1 
GeV using Eq. (D.ll) and Eq. (3.53) with different assignments of the short distance parameters. 



cuts, respectively. Given how small the effect from varying f^^ is, and its strong correlation 
with the variation of C, we do not include it in the overall uncertainty estimate for .Rcut- 

Parametric Uncertainty 

The numerical values reported above depend on two short distance parameters, m^® and 
Q^^. Our default values are m^^ = 4.19 GeV corresponding to a pole mass of = m^^ = 
4.88 GeV and af^{mz) = 0.119 corresponding to af^{mi,) = 0.217. The number of light 
flavors was set as Nf = 4. 

To estimate the uncertainty in the computed values of Rcut we repeat the calculation 
with different assignments of these parameters. The results are summarized in Table 4. 
Clearly, the largest parametric uncertainty is related to the value of the quark mass. 

Uncertainty estimates for P+ < P^a^ cut 

Let us now perform the calculation of i?cut with a cut on the small lightcone component 

in addition to a mild cut on the lepton energy (Ei > 1 GeV). As before, the central 
value is obtained by numerical integration of the differential distribution computed using 
Eq. (D.l) with the matching of Eq. (3.45) over the relevant P~ and Ei phase space, where 
the renormalization scale in the matching coefficient is set to be /x = P~ and C = 1. The 
central value computed this way for P+ < 0.66 GeV and Ei > 1 GeV (as in the Belle 
measurement [7]) is Pcut = 0.535. 

The error analysis completely parallels the one applied for the Mx cut. We therefore 
briefly summarize the results. The matching scheme uncertainty is summarized in Fig. 10. 
Numerical values for P+ < 0.66 GeV are given in Table 5. 

The renormalization scale dependence in the matching coefficient is shown in Fig. 11. 
The result for Pcut of P"^ < P^ax = 0-66 GeV (and Ei > 1 GeV) is summarized in 
Table 6. The total uncertainty due to unknown higher-order corrections amounts to ±7.8%, 
somewhat higher than the uncertainty for the Mx = 1-7 GeV cut, which is 5.6%. 

Next, the uncertainty in the shape of the quark distribution function is obtained by 
modifying C between C = 0.1 and C = 10. The result is summarized in Fig. 12 and Table 7. 
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Figure 10: The event fraction i?cut within the range P+ < P^^^ and Ei > 1 GeV, plotted 
as a function of the cut value i^niax- Three different matching schemes are used: (a) the partially 
integrated distribution with respect to £"/, namely Eqs. (D.8) and (D.9) with Eq. (3.53), numerically 
integrated over P~ (full line); (b) the fully differential distribution computed by Eq. (D.l) and the 
matching of Eq. (3.45), numerically integrated over the Ei and (dashed line); and (c) same as 
in (b) but with the matching of Eq. (3.40) (dotdashed line). 



matching scheme 


Event Fraction Rent (-P+ < 0.66 GeV, Ei > 1 GeV) 


variation (%) 


(a) 


0.5118 


-4.4 


(b) 


0.5353 


default 


(c) 


0.5659 


+5.7 



Table 5: Results of the event fraction i?cut within the range P+ < P^^^x — 0-66 GeV and Ei > I 
GeV using different matching schemes (see Fig. 10 and discussion next to Table 1) with the scale 
of the coupling in the matching coefhcient set as fj, — P^ . 



Finally, the parametric uncertainty estimates are summarized in Table 8. The total 
parametric uncertainty is ±12.5. In conclusion, when applying a P"*" < 0.66 cut instead 
of a P+P- < Ml- = (1.7GeV)2 cut the event fraction is reduced by 13% - 14%, while 
the theoretical uncertainty increases. This applies separately to all sources of uncertainty, 
namely the unknown higher-order corrections (NNLO), the shape of the quark distribution 
function (C) and the values of as and rrn,. 

4.3 Extracting \Vuh\ 

Let us consider first a cut on the invariant mass of the hadronic system, as well as a mild 
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Figure 11: The event fraction _Rcut within the range P+ < -P,iax Ei > 1 GeV, computed 
using Eqs. (D.8) and (D.9) and Eq. (3.53), plotted as a function of the cut value Pj^ax- Different 
assignments for the argument of the coupling in the matching coefficient are used. 



scale of Os 


Event Fraction Rent for P+ < 0.66 GeV and Ei > 1 


variation(%) 


f, = P~ 
fj, = mb 
fj, = mb/2 


0.5118 
0.5385 
0.4847 


default 
+5.2 
-5.3 



Table 6: Results of the event fraction within the range P+ < 0.66 GeV and Ei > 1 GeV using 
Eq. (D.9) and Eq. (3.53) with different assignments for the scale of the coupling in the matching 
coefficient (the n — 1 moment). 



c 


Event Fraction Rent {P^ < 0.66 GeV, Ei > I GeV) 


variation (%) 


0.1 


0.4971 


-2.9 


1 


0.5118 


default 


10 


0.4909 


-4.1 



Table 7: Results of the event fraction within the range P+ < PjJax = 0-66 GeV and Ei > 1 GeV 
using Eq. (D.9) and Eq. (3.53) with different assignments of C in Eq. (3.29). 

cut on the charged lepton energy, Ei > Eq, as in Ref. [7]. The calculation of the partial 
branching fraction takes the form: 

AB{B — > XJv,P+p- <Ml,Ei> Eq) 
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Figure 12: The event fraction i?cut within the range P+ < -P,uax Ei > 1 GeV, computed 
using Eqs. (D.8) and (D.9) and Eq. (3.53), plotted as a function of the cut value P,^ax- Different 
assignments for the parameter C, which controls the shape of the perturbative quark distribution 
function (see Eq. (3.29)) are used. 



parameter 


Event Fraction R^ut {P^ < 0.66 GeV, Ei > 1 GeV) 


variation (%) 


af^{mz) = 0.115 


0.4922 


-3.8 


af^{mz) = 0.119 


0.5118 


default 


af^{mz) = 0.122 


0.5209 


+1.8 


^/ = 3 


0.5207 


1.7 


Nf = 4 


0.5118 


default 


= 4.14 


0.4513 


-11.8 


mf^ = 4.19 


0.5118 


default 


mf^ = 4.24 


0.5650 


+10.4 



Table 8: Results of the event fraction within the range P'^ < 0.66 GeV and Ei > 1 GeV using 
Eq. (D.9) and Eq. (3.53) with different assignments of the short distance parameters. 



= tbT,^,^, {B XJi?) RcntiP^P- < Ml, El > Eo). (4.8) 
The total effect of the cuts amounts to the following event fraction: 

r(P+p- < (1.7GeV)2, El > lGeV)/ro 
T{P+P- < Ml, El > 0)/ro 

0.615 ± 0.035[higher orders] ± 0.014[quark distribution] ± O-^^^ [parametric (mfS)] ' ^'^'^^ 



cut 
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Here the central value was computed using numerical integration of the fully differential 
distribution, Eq. (D.2) with the matching of Eq. (3.45) where the renormalization scale 
was set to /i = P~. The theoretical error estimate was done as explained in Sec. 4.2: 
the uncertainty associated with higher-order corrections is ±5.6%, the one associated with 
the quark distribution function (renormalons on and power correction on the soft scale) 
is ±2.4% and the parametric uncertainty (dominated by the uncertainty in rrib) is ±7.4%. 

Using Eq. (4.9) together with the theoretical result for the total charmless semileptonic 
width of Eq. (2.1) in Eq. (4.8) with the experimental value of the B lifetime tb = (1.604 ± 
0.011) ps, and comparing the result to the Belle measurement^^ [7], 

AB{B — ^ XJi?,P+p- < (1.7GeV)^ Ei > IGeV) = 1.24 • lO'^ (±13.4%) (4.10) 

we get, 

= (4.35 ± 0.28[e.p] ± 0.14j^^_^^^^^(^p^, ± 0.22[,,_,„,,j) • lO'^, (4.11) 

where the three sources of errors quoted separately are: (1) the total experimental error [7] 
on the measured BF in Eq. (4.10); (2) the error on the total width, dominated by the 
uncertainty in in Eq. (2.11); (3) the theory error on the event fraction associated 
with the hadronic invariant mass cut, computed by adding the three sources of uncertainty 
in Eq. (4.9) in quadrature. 

Let us consider now the cut on the P"*" momentum as proposed in [16] and measured 
in Ref. [7] with P+ < = 0.66 GeV. The calculation of the partial branching fraction 
takes the form: 

AB{B XJu,P+ < P+,,, El > Eo) 

= TsAotai {B Xulv) Pcut(P+ < P+ax, El > Eo). (4.12) 

The effect of the cuts amounts to the following event fraction: 

_ r(P+ < 0.66 GeV ,Ei>l GeV)/ro 
" r(P+ < Mb, El > o)/ro 

= 0.535 ± 0.042 [higher orders] =•= 0-021 [quark distribution] =t ^■^^^[■pg^rametric {mf^)Y^^-^^^ 

Comparing Eq. (4.12) to the Belle measurement [7], 

AB{B — > XJv, P+ < 0.66 GeV ,Ei>\ GeV) = 1.10 • 10"^ (±17.2%), (4.14) 

we get, 

iKbl = (4.39 ± 0.36[e.p] ± 0.14j^^_^^^^,(^Ms^j ± 0.38[,h-cuts]) • 10"^ (4.15) 

Comparing \Vuh\ from the P"*" < 0.66 GeV cut in (4.15) to the one from P^P~ < 
M\ = (1.7GeV)2 in (4.11) and with Ref. [7] we find: 

'^^The experimental uncertainty is computed by adding the systematic and statistical errors in quadrature. 
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• The central values of we extract using the two cuts agree very well. 

• Our results also agree very well with the value quoted by Belle [7] based on their 
own analysis of the Mx cut, while the central value obtained there for the P"*" cut 
is ~ 5% higher^'^. 

• The P^" cut is characterized by larger theoretical uncertainty; this is even more pro- 
nounced than the increase in the experimental error. The largest sources of theoretical 
uncertainty, for both cuts, are the value of the quark mass and the unknown NNLO 
corrections. 

5. Conclusions 

We have presented a calculation of the fully differential distribution in B — > X^W decays 
that relies only on resummed perturbation theory, and yet describes the meson decay 
spectrum with approximately correct physical support in the small -P+ region. What makes 
this at all feasible is that the leading renormalon divergence of the Sudakov exponent is fully 
under control; the corresponding ambiguity of the resummed moments, which amounts 
to a global 0{A/P~) shift, cancels out exactly upon converting to hadronic variables, 
a conversion that itself involves the a priori ambiguous pole-mass definition. While a 
regularization (Principal Value) has been arbitrarily chosen at the intermediate stage, the 
end result is unambiguous. 

In sharp contrast, an attempt to use Sudakov resummation with fixed logarithmic 
accuracy — no matter how high — is bound to fail dramatically due to the unregulated 
renormalon divergence^ ^ . An example is presented in Figure 13 that compares the DGE 
spectrum at Ei = 1.5 GeV and P~ = 3.5 GeV to that obtained by standard Sudakov 
resummation with LL, NLL and NNLL accuracy. Obviously, DGE makes a qualitative 
difference. 

It should be stressed that the detailed description of the spectrum in the P+-peak 
region, and especially near the P^ = endpoint, is necessarily approximate and it must 
rely on certain assumptions. First of all, there are both perturbative and non-perturbative 
contributions. Using moments the latter appear as powers of 0{nA/P~), so they are 
parametrically small at small n, but they gradually become important as n increases. 
Uncertainties corresponding to the first and hopefully the second powers of nA/P~ can 
be avoided — but not the high powers. Towards the small-P+ endpoint {P~^ — > 0) high 
moments become relevant and the parametric hierarchy is lost. What dictates the details 
of the decay spectrum in this region is the quark distribution in the meson. Physically 
one expects that this distribution would differ from the one in an on-shell quark owing to 
the interaction between the b quark and the light degrees of freedom. It is our working 
assumption that these differences are small and therefore better calculation of the latter 

^"The values corresponding to the two cuts there are still consistent within errors. 

•^^ Moreover, in this case it is not at all obvious which mass should be used in the conversion of the partonic 
spectrum to hadronic variables. 
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0.20 




Figure 13: The differential (left) and integrated (right) spectrum based on DGE and on the 
fixed-logarithmic-accuracy formula of Eq. (3.26), plotted as a function of P+. All calculations 
are matched to NLO. The LL, NLL and NNLL accuracy results are plotted as dotdashed, dashed 
and full line, respectively. The DGE curve smoothly extends below the partonic phase-space limit, 
filling up the region < < A. In contrast, the three fixed-logarithmic-accuracy curves end 
above A, at the point where the resummed results become complex owing to the Landau singularity. 



amounts to better description of the former. Eventually the difference between the two 
should be quantified. 

Non-perturbative contributions aside, computing the on-shell decay spectrum itself 
within a given regularization scheme for the renormalons is still an incredible task. In 
Ref. [31] and in the present work we made significant progress in this direction. As we saw, 
the main difficulty is that the Borel sum in Eq. (3.43) is not entirely dominated by the 
perturbative expansion in powers of u; additional information on the behavior of the Borel 
integrand away from the origin is needed. Since the general structure of the exponent as 
a function of n is known, the main uncertainty arises from the soft anomalous dimension 
function Bs{u). The cancellation of the leading renormalon ambiguity with the pole mass 
allows us to determine Bs{u = 1/2) accurately [31]. Based on the known analytic result 
in the large-/3o hmit, we further assumed here that Bs{u = 1) = and that this function 
has no other zeros. We then parametrized the contribution from the region Bs{u ^ 3/2) 
to facilitate uncertainty estimates. We found that approximately correct non-perturbative 
support properties are recovered in the Principal Value prescription if this contribution is 
indeed small. This is demonstrated in Fig. 8. 

Under these assumptions we obtained definite predictions for the B — > [31] 
and B — > Xylv spectra. So far comparison with experimental data has been performed 
for the first two central moments in B — > Xg^ with varying cuts, finding very good 
agreement [41]. More detailed comparison of both distributions has a great potential in 
testing these assumptions and, eventually, in quantifying power corrections. 

The fully differential B — > Xylv spectrum computed here can be used to extract 
\Vub\ from any experimental measurement of the partial branching fraction. The resum- 
mation employed is specifically designed to improve the determination of the partial rate 
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in the small hadronic mass (or small P^) region, which is the experimentally-favorable 
way of discriminating charm. By numerically integrating the differential spectrum over the 
relevant phase space we computed the fraction of events measured for a given experimen- 
tal cut. Using recent measurements by Belle [7] for the branching fraction in the regions 
P+P- <M\ = (1.7GeV)2 and P+ < P^^ = 0.66 GeV (with Ei > 1 GeV), we extracted 
\Vub\ and analyzed the theoretical micertainty. The final results are quoted in Eqs. (4.11) 
and (4.15), respectively. We find that 

• The most important source of uncertainty is the value of the b quark mass. It 
directly influences both the calculation of the total rate and that of the event fraction 
associated with the cut. Better determination of the mass will therefore translate 
directly into better measurement of \Vub\- One natural possibility for the b quark 
mass measurement is to make use of recent data on the B — > Xgj spectrum [12, 13] 
together with calculated spectrum [31]. 

• Another important source of uncertainty is higher-order perturbative corrections. 
NNLO calculation of the fully differential B — > XylU is not beyond reach. When 
available, it can be readily matched to the resummed spectrum. This will signifi- 
cantly reduce the matching scheme dependence as well as the renormalization scale 
dependence in the matching coefficient. 

• Finally, under the assumptions discussed above, the uncertainty associated with the 
detailed shape of the quark distribution function is found to be quite small. This is 
most striking in the Mx cut, see Fig. 9 and Table 3. 
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A. Resummation in terms of hadronic and leptonic invariant masses 

Let us choose the following kinematic variables [30] , all defined at the partonic (perturba- 
tive) level: 



where momentum conservation implies ph = Pj +q, where pj is the total momentum of the 
jet and q = ki + is the total momentum of the leptons. The total semileptonic width 




(A.l) 
(A.2) 
(A.3) 
(A.4) 
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into u-quark is: 



ftotai = / dz I dx 







(ir(z,^,x) 



1 + 



CFas(m) / 25 vr 



TT 



/a;-2(i-i/x) dzdxdC 

where Tq is given by Eq. (3.6). 
NLO result 

The perturbative expansion of the triple differential width take the form: 

1 dT{^,z,x) 



+ 



(A.5) 



To d^dzdx 



V{z,x)5{0+n{i,z,x) 



where 



V{z, x) = vo{z, x) + ^^'^^^^^\ i{z, x) -\ , 



TT 



TZ{(,,z,x) = ^^^^^^—^di{^,z,x)-\ . 

vr 

The NLO result for the virtual coefficients is: 
vo{z, x) = 12{z + 1 — x) {x — z) 

vi{z,x) = -12(^-2In(l-z)2-^7r2-^-Li2(z)^ {x - z - l){x - z) 

— 61n(l — z){^x — hz — A){x — z) 

while the real-emission coefficient is: 

' Q{x - z - C) {z - x)t{z + i - 2) 



(A.6) 

(A.7) 
(A.8) 

(A.9) 



di{i,z,x) = 



((z + 02-4z)2(l-0 
{2z^ -lOz^ + 4:z'^C-'^0z^ + 3^^ z + 13z + f - 30 

6 2 < + C - 2) (5 - G r + 2 r ( + + 2 5^ - 1 



(A.IO) 



+ 



iz + ^-2) 



,, ^ ,„ ^{z - xz - x^^- X )t 

{{z+if-^zf 

2z'' + 6z^^- 10z^ + ^2 + 21^2 _ 28^2^- 14^2^- lOz 



+ 4zC^ + 24z,^-5^^ + 2^3^,t^ 



1 



In 



f + t 



3(7^2 +4z,^- 18^ + 6^ + ,^2) (^_^_^) [-z^x) 

(-1 + ((^ + 0^-4^) 

3 z (-9 z + z ^ + 10 + ^2 - 3 3{z-xz-xi^x^) 

^ ((z + 02-4z) + ((^ + ^)2_4^)2 ^ 

(^1 1 + 17 ^ - 28 + 7 ^2 ^ + 20 z - 38 z C + 10 ^2 + ^3 j 

with t = \J{z + 0^ ~ 4^;/(2 — z — 0- The ^ — >■ 1 singular (non integrable) terms in di, 
namely 

3(1 + z - x)(x - z) _ ^ g j^^^ J 



1-e 



are regularized as ()* distributions, as in Eq. (3.11). 
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Resummat ion 



The all-order structure of the ^ — > 1 singular terms in the triple differential width in the 
large-/3o limit has been determined in Sec. 3 of Ref. [30] using the Borcl technique. The 
result is summarized by Eq. (3.15) above. Sudakov resummation can now be performed by 
exponentiating these singular terms in moment space. Defining moments by 



dTN(z, x) 
dzdx 



1 



r 



Jx-z(l-l/x) 



dec 



N-idT{^,z,x) 



d^dzdx 



(A.12) 



To Jo 



N-i dT{^,z,x) 



d^dzdx 



+ 0{1/N) 



The singular terms for ^ — > 1 generate terms that contain powers of In A^, while 0(1/N) 
terms are neglected. After matching the resummed distribution can terms can be obtained 
by the following inverse Mellin formula: 



1 dT{^,z,x) 



'^+^°° dN drNiz,x) 



To dzdxd£, Jc-ioo 27rz dzdx 
The matched moment-space result takes the form: 



-N 



matched 



dT]\f{z, x) 



exp 



dzdx 
Cf 



matched 

°^ du 



V{\x)+KD,\x,x)^^^^ + 0{al) 



TT 



(A.13) 



(A.14) 



/3o Jo u 



T{u) 



A' 



2 \ " 



Bs{u){l-zf''V{-2u) 



T{N) 



1 



r(A - 2u) r(i - 2u) 



To perform matching at 0{as) one needs to compute the moments at this order 



<(A,^) ^ ^ d^e 

-L Jx-z{l-l/x) 



d^dzdx 



di{^,z,x), 



where di{^,z,x) is given by Eq. (A. 10). We find: 

1 



(A,x) =vo{z,x) 



*i(A)--7r^-(7E + ^(Ar))^ 

D 



(A.15) 



(A.16) 



+ 



7 



21n(l-z) {^iN) + jE) } + 0(1/Ar). 



These terms, including the constant for N — > oo, are reproduced by the expansion of the 



-N 



exponent in Eq. (A.14), so ADi (A, x) contributes only at the level of 0{1/N) corrections 

- — ^ — ^oo 

{ADi (A, x) = 0). Therefore, Eq. (A.14) with the coefficients of V(A, x) given by Eq. (A. 9) 
can readily be expanded to predict the log-enhanced terms at 0{a^). Converting this result 
back to ^ space, we get: 



1 dT{^^,z,x) 



Tq dzdxd^ 



vo{z,x) + 



In(l-g) 
1-e 



+ (-^ + 21n(l-z; 



1 
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vo{z,x) + vi{z,x) 



+ 



TT 



Cl 

2 V 1-e 



(A.17) 



+ 



21 



31n(l-z) + -JC^ + 
^2 49 



2 , ^ll^A_iV^)^^ 



+ ( (4 1n(l-zf -^ + ^-71n(l-z))c2 + ( 



13 
'72 



+ 



Nf + 



6 16 



2 
3 

ln(l-0 



1 



ln(l - z)7r^ 



1 



32 



, 17 , , 177r2 1 11 , ^2 

+ — ln(l - z)^ h - Ca H ln(l - zY 

^9 ^ ^ 72 4 6 ^ ' 288 



i ln(l-z)7r2 ) CfCa 



1 



1-e 



+ 



+ 2 ln(l - 2) 



1-e 



vo{z,x) + 



/In(l-e) 



1-e 



TT ) 



Changing variable in Eq. (A.17) to the hghtcone momenta p and A we get 



1 cir(A,/9, x) 



To dXdpdx 



= vo{z = 1 — A, x) + 



lnp\ 



+ ln(A) - 



7 



vo{z = 1 — A, a;) + vi{z = 1 — X, x) 



P J 



TT 



+ 



21 Ci 



+ [ --Cl ln(A) + 
+ ( -ic|ln(A)2+ ( - 



+ 



7Cl 



+ 



Nf 11 Ca\ 

6 12 y 



P 



Cf ln(A) 



(A.18) 



+ (^C|ln(Af + 



35 Cl 



+ 



CpNf 



"6 



2 1 

2 



3 

32 



Cs - TTTT ) + 



36 144 



1 177r2 905\ 
+ (-C3 + ^-^j CfCa 



vo{z 



X,x) 



+ 



Inp 



+ ln(A) 



7 



= 1 - A,x; 



TT 
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These logarithms agree exactly with those predicted using the variables r and A, Eq. (3.5) 
with Eq. (3.41). Note, however, that the separation into virtual and real is different when 
different variables are used to define the ()* distributions. 



B. Calculation of the real-emission contribution at NLO in moment space 



Singular part 

For the singular part of Eq. (3.32) we get 



K[''"'<'{\,x) = -woiX,x) dril-r) 

Jo 



n—l 



ln(r) + | 



(B.l) 



= -wo{\,x) 



where we applied the definition of the ()* distribution in Eq. (3.11) and defined 



ioifm) = / dr 



(1 - r)"-^ - 1 



Inr 

r 



(B.2) 



= (n - l)r„ 



4F3([1, 1, 1, 2 - n], [2, 2, 2], r^) - ln(r„) 3F2([1, 1, 2 - n], [2, 2], r„ 



TT 



-\n{rm)si{n, 1 - r^)-S2(n, 1 - Vm) 



and 



Jo 



dr 



J r 

r^(n-l) X 3F2([l,l,2-n],[2,2],r„) 
In(r^) + *(n) + 7E - si(n, 1 - r„,) , 



(B.3) 



where si^2 denote the infinite sums: 

~" (l-r^)"+'= 



si(n, 1 - rm) = 



k=0 

(1 - rmT 



n + k 



(B.4) 



n 



2Fi([l,n],[n + l],l-r„) = (1 - r„)" LerchPhi(l - 1, n) 



S2 



(n, 1 - r„) = ^ 



fe=o 



(1 - rm) 
n + k 



n+k 



*(n + k)- *(n) 



r„)«r(n) 



Alternative expressions, useful for small Vm, are: 

oo 

Slin, l-rm) = -{ Hrm) + *(n) + 7.) - g , ^ (n - k) 

S2{n,l-rm) = l{Hrm) + nn)+lEy + ^^'-\^i{n) +EmTT"^'^^"^ 



(B.5) 



^^r{k + l)k^r{n-k)' 
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iFvom the definitions of io,i and Eq. (B.l) it follows that the first moment {n = 1) of 
the singular part, which is relevant for the total width, is simply: 



K^-^'"'^{X,x) = -wo{X,x) 



1 7 
-ln2(r„) + -ln(r„) 



(B.6) 



We note that ^qi ~ C>{rm), so for any moment the small-^^ limit {x — > 1) is 
particularly simple: 

ifi"^'°« (A,x ^ 1) Xp^^'°^ (A,x 1). 

For fixed the sums si^2('^) 1 — fm) in Eq. (B.4) vanish exponentially at large n. Thus, in 
the final expression in Eqs. (B.2) and (B.3) only the square brackets contains terms that 
do not vanish at n — > oo. The logarithms of appearing explicitly in Eq. (B.l) cancel 
against the ones in and io, yielding 



K'l-''^iX,x) = -wo{X,x] 



^-^'i(n) + (^'(n)+7Ef 

D 

7 



--{^{n)+^E) (B.7) 



-si{n, 1 - rm) ln(r„i) - S2{n, 1 - r^) - -si(n, 1 - r„ 



so the r^n dependence of K""'"^ (A,a;) appears only though 0{l/n) terms. The terms is the 
curly brackets in Eq. (B.7) in the first and second lines can interpreted as originating in 
the two following integrals, respectively, 



Kr--iX,x) = -woiX,x)!^ J\l - vT-' + 



(B.8) 



•Jr-m. 



1 -r) 



1 / In r 71 

- + 7- 
r 4 r 



Regular part 

Let us consider now the regular part of Eq. (3.33). The first moment (n = 1), which is 
useful in calculating the integrated width, is fairly simple: 



K^-^ "''■{X,x) = 6Li2 ( 1- 



1-x 
X 



8 A + 36 + 232 A'^ - 190 A^ - 90 A^ 



(B.9) 



- 4 - 196 A^ - 76 A + 225 A^ + 40 aM (1 - x) 



+ (34A - 33 A^ + 8A^ - 8) (1 - xf 



In 



1 — a;\ 1 — X 



216 A^ - 1140 A'' 



+ 84A + 1416A2 -612A2 -3 - 199 A^ + 211 A^ + 2 - 70 A + 44 A"^ (1 - a;) 



2(1-A) A2-2A + 2 (1-x 



+ 



6A 



540 vr^ A^ - 216 tt^ A^ - 1392 tt^ A^ 



+ 11407r2 A5- 487r^ A 



2 \2 



+ {l-x) 



4 TT^ - 516 A - 196 TT^ A^ - 76 tt^ A 



-52- 



+ 225 TT^ + 216 A"^ - 1140 A^ + 1365 A^ + 60 + 40 tt^ A^ 



{l-xf 



-624 A 



2A 

+ 39 + 16 TT^ A"^ + 1995 A^ - 1812 A^ + 372 A"^ - 16 tt^ A - 66 tt^ A^ + 68 tt^ A^ 



29A2 -70A + 43 



To express a general n-th moment of the regular part, we write 

2 2 



kr'\\r,x) = /i(A,x)(l-r)Mn(r)+ J] /,(A, x)(l - r^, (B.IO) 

i=-4 'i=-3 

where we consider n complex with sufficiently large real part, Re(n) > 3, to obtain conver- 
gence of the defining moment integral in Eq. (3.31) for individual terms in the sum, and 
where 

/_4(A, x) = -36 (A - 6) (A - if (A - 1 + xf (B.ll) 
/_3(A, x) = 36 (A - 1) (3 A^ - 17 A + 9) x^ + 36 (7 A^ - 40 A + 18) (A - 1)^ x 

+ 36(A-1)^(4A2-23A + 9) 
/_2(A, x) = -6 (19 A^ - 102 A^ - 34 + 112 A) x^ 

- 12 (-1 + A) (29 A^ - 159 A^ + 148 A - 38) x - 4008 A^ - 1608 A 
+1812 A^ - 240 A^ + 252 + 3792 A^ 

/_i(A, x) = 6 (8 A^ - 33 A^ - 8 + 34 A) + 6 (-212 A^ + 291 A^ - 144 A + 20 + 40 A^) x 
+ 2616 A^ + 708 A - 1380 A^ + 216 A^ - 72 - 2088 A^ 
/o(A, x) = -6 (A - 1) (a2 + 2 - 2 A) x^ - 6 (-58 A^ - 34 A + 71 A^ + 6 + 14 A^) x 

- 990 A-^ + 24 - 114 A^ + 600 A'' - 228 A + 708 A^ 

/i(A,x) = 12(A-1)A(A2 + 2-2A)x + 6(A-1)A(6A^-19A2 + 18A-6) 
/2(A,x)= -6A2(A-1)(A2 + 2-2A) 



and 



/_3(A, x) = -36 (A - 6) (A - 1)2 - 72 (A - 6) (A - 1)^ x + 216 + 360 A^ 

- 36 A^ + 1440 A^ - 900 A - 1080 A^ (B.12) 
/_2(A, x) = 18 (A - 1) (5 A^ - 27 A + 12) x^ + 108 (2 A^ - 11 A + 4) (A - 1)^ x 

+ 126 A^ - 1080 A"^ + 2700 A^ - 216 - 2880 A^ + 1350 A 
/_i(A, x) = -6 (12 A^ - 58 A^ + 60 A - 19) x^ 

- 6 (A - 1) (41 A^ - 209 A^ + 174 A - 46) x + 162 - 1008 A + 1302 A^ 
+ 2472 A^ - 2748 A''^ - 180 A^ 

/o(A,x) = 3(A-1)(6A2-14A + 7)x2 

+ 3 (21 - 120 A - 210 A^ + 262 A^ + 44 A"^) x - 42 + 342 A - 804 A'' 

- 1053 A^ + 1422 A^ + 135 A^ 
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/i(A,x) = -6A(5A-7)(A-1)2x-63A + 261A^ + 252A2-396A3-54A^ 
/2(A,x) = 3A2(3A-7)(A-1)2 

Despite appearance ki'"'^'{X,x) is regular also for r — > 1 (recall that the phase space 
is r < (1 — x)/A < 1): there are cancellations between the would-be singular terms in 
Eq. (B.IO). 

Performing the integration in Eq. (3.33) we obtain the following moments: 



K^"'^ {X,x) = ^ fi{X,x)In+i{rm) + 



(B.13) 



i=-4 



i=-3 



where we define: 



and 



Jn{rm) = / dr{l - r)"-i = - (1 - (1 - r„)") , 
Jo 



11 

In{rm) = / dr{l - r)"-i Inr = ^- -{^{n) + 7e) + 



(1 - rmT+^ 



- (1 - (1 - r^r) + (1 - rmT Hrm) + (*(n) + 7e) + £ " '''"^ 



(B.14) 
(B.15) 

n+fe " 



n 



k=o 



n + k 



i^(r^) + (l-(l-r^n{ln(r^)-i 



where ii{rm) is given in Eq. (B.3). Note that if" (A, a;) has no singularities^^ on the real 
n axis, neither at positive nor negative values of n. This follows directly form Eq. (3.33) 
as the upper limit of integration obeys r„i < 1. Note also that for fixed n the functions 
Jn{fm) and Inifm) Vanish in the r„j — > limit, and the leading order approximation to 
these functions in this limit is n-independent: 

Jn{rm) ^ Jn=l{rm) + C(''m); Jn=l{rm) = Tm (B.16) 

IniTm) - In=i{rm) + 0{r'^)] In=i{rm) = r^(ln(r^) - 1). 



C. Matching coefficients for the electron— energy integrated width 

Here we compute the matching coefficients Vr{X, x > xq) and AKi\X, x > xq) in Eq. (3.53), 
by integrating the virtual coefficients of Eqs. (3.7) and (3.8) and the real-emission moments 
of Eq. (3.46), respectively. 
The virtual terms are: 

Vr{X,X > Xo) = Wo{X,Xo) + ^^"^^'"^^ Wl{X,Xo) H , (C.l) 

TT 

^^In Eq. (B.15) there is exact cancellation of singularities at all non-posituve integer between the ^ term 
and the sum. 
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-r(A) 


12A2(1 - A) 




12A2(-1 + 2A) 




-12A3 



Table 9: The Born level coefficient, wo{X,x) of Eq. (3.8), decomposed according to Eq. (C.5). 



where the LO coefficient is 

'w^{X)= 2A2(3-2A), A<l-a;o 



w^o{\,xo) = 4x(] + 6 (-3 + 2A) + 12 (1 - A) (2 - A) xo 



-10-12A^ + 24A, 



X>1 — Xq 



and the NLO coefficient is: 



wi{X,xo) 



w\{X,xo) 



(C.3) 



with 



X < 1 — Xq 
X> 1 — Xq 

<(A) = 2A2(-3 + 2A)Li2(l - A) + A^(4A - 9) ln(A) + ^X^{47r^ + 15)(-3 + 2A) 

w'i{X,xo) = 2(l-a;o)(2xg-7xo + 6xoA-12A + 6A2 + 5)Li2(l-A) (C.4) 
+(1 - xo)(-17xo + ixl - 30A + 12A^ + 13 + 12xoA) ln(A) 

+^(1 - a;o)(47r2 + 15)(2x^ - 7xo + 6xoA - 12A + GA^ + 5) 

Eq. (C.2) shows that the small A limit is power suppressed. This property is crucial for the 
validity of the perturbative treatment, as the effective hard scale in the Sudakov factor is 
Xmi,. 

The integration over AK-^ (A, x) is more involved as it depends on x in a non-trivial way 
through = (1 — x)/X, which was the upper limit in the r integration. It is convenient 
to arrange the polynomial dependence on x in powers of rm, namely 



«;o(A,x) = 5;4^\A)r^; /,(A, x) = J] /(''^^(A) ; /,(A,a;) ^ J] /(^'^^(A) r^, (C.5) 

j=0 j=0 j=0 

where the coefficients are given in Tables 9 and 10, and change variables from x to 

i; = l-rm = l-(l- x)/X. 

One then obtains: 

AKr{X,x>xo)= [ dxKKl{\x) = xj dvKxliX, x = I - X{1 - v)) (C.6) 
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i 


J 


/(*,j-)(A) 


/(i,i)(A) 


-4 
-4 
-4 



1 

2 


— 36(A — 6) (1 - Aj^A'' 
72(A - 6) (1 - A)2a3 
-36(A - 6) (1 - A)2a3 







CO CO CO 

1 t 1 



1 
2 


-36A''(1 — A)(4A^ — 24A + 15) 

36A3(1-A)(7A2-41A + 24) 
-36A3(1 - A)(3A2 - 17A + 9) 


— 36(A — 6)(1 - Aj^A"* 

72(A-6)(1-A)2a3 
-36(A-6)(1-A)2a3 


-2 
-2 

-2 



1 
2 


-6X"^{8 - 120A + 311A'^ - 244A'^ + 40A^) 
12A2(4 + 205A2 - 74A - 169\^ + 29A^) 
-6A3(-34 + 19A=^ + 112A - 102A2) 


-18A"*(1 - A)(7A'^ - 41A + 24) 

72A3(1 - A)(3A2 - 17A + 9) 
-18A3(1 - A)(5A2 - 27A + 12) 


-1 
-1 
-1 



1 
2 


1 c\ \ "2 / A A f \ 1 -i-i/^\2 rir'X.S 1 -io\4\ 

12A (4 — 45A + 116A — QSA"* + 18A ) 

-6A2(4 + 225A2 - 76A - 196X^ + 40A^) 
6A''^(-8 + SA''^ + 34A - 33A2) 


— 6A (8 — 87A + 220A — 176A'' + 30A ) 

6A2(8 + 267A2 - lOOA - 226A3 + 4^"^) 
-6A3(-58A2 + 60A - 19 + 12A3) 








1 
2 


— 6A^(8 - 50A + lOBA^ — 86A-' + 19A^) 
6A2(2 + 65A2 - 26A - 56X^ + 14A^) 
-6A=^(-1 + A)(A2 -2A + 2) 


3A^(15 — 109A + 270A — 224A^ + 45A^) 
-3A2(7 + 222A2 - 78A - 198A3 + 44A^) 
3A3(-1 + A)(6A2 - 14A + 7) 


1 
1 
1 


(J 

1 
2 


— oA (1 — Aj(DA — 17A + 14A — 2) 
12A3(1 - A)(A2 -2A + 2) 



0\2/1 W/ 1 1 O \ \ /A \ 2 \ 1 T\ 

oA (1 — Aj(— I + zAj(9A — zoA + 7j 
6A3(l-A)2(5A-7) 



2 
2 
2 



1 
2 


6A3(1 -A)(A2 -2A + 2) 




3A3(3A-7)(1 -A)2 





Table 10: The NLO matching coefficients, MX,x) and /i(A,a;) of Eq. (B.ll) and Eq. (B.12), 
respectively, decomposed according to Eq. (C.5). 



where vq = fo(A, xq) = I — (1— a;o)/A. The integration in Eq. (C.6) can be done analytically. 
Let us write the result as 



AKi (A, X> Xq) = < 



^Ci^(n,A,i;o(A,a;o)) =Cp' + q*^ 



A < 1 — xo 
A > 1 — a;o 



The expressions for the first moment (n = 1) are simple (see Eq. (3.47)): 

Cp-^(n = 1, A,i;o) = (-3A2 + 3A\g - IX^v^ + 2A^) ln2(l - vq) 



(C.7) 



15 



17, 



+ ( (4A^ - 6X')vo +(2X' + ^X' ) - ^X'v^ - ^X' - ^X' ) ln(l - ^;o) 



17 



+^A\3 + ( ^A-^ _ z_x^ + ( __A^ _ _A^ ) + -_x^ _ r_A^ , 



9 



11 



15 



6 



1 



21 
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18 



\n = l, X, vo) = -X^ - 16 A - 66 A^ + 16 A^ + 68 A^ - 180 A vq 

+72 A^ Vo - 390 A^ vq + 12 + 471 A^ vo + 24 A'' + 246 A^ - 132 A + 24 
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-162 )U2{vo) + ^A^ (1 - vo) yl2v^ - 9v^ + 3v^ - 6v^ X 

v4„,2 ionn\3„,2 , io/?«\2„,2 aoa \ „.2 , io„,2 i ott \4 



+287 A* - 1299 A'* + 1366 A^ - 434 Xv^ + 12v^ + 377 A* 
+3301 A^ vo + 264 vq - 2361 A^ - 1586 A + 5 A^ + 31 A^ + 58 A - 39 A^ 

-60 ) ln(l - vo) + ^A^ (A - 1) (61 A^ - 146 A + 92) + ^A^ (1229 A^ 
J 24 18 y 

-6141 A^ + 6673 A^ - 2048 A + 129 + 7r^(48 A^ - 198 A^ + 204 A^ - 48 A)^ 
+^A2 f - 385 A^ + 144 TT^ A^ + 1077 A^ - 780 tt^ A^ + 199 A^ + 942 tt^ A^ 
- 1 196 A - 360 TT^ A + 390 + 24 TT^^ + ^ A^ 1^ - 233 A** + 24 vr^ A^ + 1563 A^ 
- 162 TT^ A^ - 2400 A^ + 246 tt^ A^ - 132 tt^ A + 1368 A - 303 + 24 tt^ 

+^ l72^ 24^ l8"^ + "6"j' 



and 



while for general complex n (with Re(n) > 3 i.e. Re(n + i) > —1, so that convergence of 
the moment integral for individual terms in the sum is guaranteed) C\{n, X,vo{X,xo)) can 
be written in terms of a few known integrals: 



C\{n,X,voiX,xo)) = J2 



2 r ^ 7 \ • 

KnjM + ^Ln,j{vo) + Mn,j{vo)] X W^'\X) 



j=0 



2 T / . 

+ E /^'''HA)^ -{nn + i + l) + lE)Pj{vo) + ^Qn+i,jM (C.9) 

-Rn+i,j{vo) - Ln+ij{vo)^ + E /^''^^(-^);^ (^Pjivo) - Qn+ijivo)^ 



All these integrals can be expressed as power series in vq, which is convenient for numerical 
evaluation: 

P,{vo) ^ f dv{l - vy = (C.IO) 

/•I d 

Rn,jM = / dv{l-vyv''\n{l-v) = —Qn,j{yo), 
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„1 oo 



oo h 



T{n)T{j + i: 



fc=0 
fe=0 



n 



j + 



+ fe (j + l)r(n + i + l) 

r(n)r(fc + l+j + n) 1_ 

r(j + 1 + n)r(A: + n + 1) n + kj' 

d 



/•I °° A; J 

K„,,(^;o) = / dt;(l-.;)Vln(l-^)^-^ = 

„1 oo 

Mn,j{vo)= / dt;(l-v)VV^— (*(n + A;)-*(n)) 

^^0 fc=o + 



r(l + n + j)(l +j)2 1^ n + A; ^ ^ ^ ^ ' 

These series converge within the physical region, i.e. < vq < 1, and complex n with 
Re(n) > 3 (which can have large imaginary part). The series expansions of Rnj{vo) and 
Knj{vo) can be readily obtained by differentiating the expressions for Qnjivo) and Lnjivo), 
respectively, with respect to j. For integer values of j, Mn,j{yo) becomes a single infinite 
sum since Qn,j{vo) reduces then to simple algebraic expressions: 

^".0 = (C.ll) 

1 V 



Qn,l = V 



1+n 



Qn,2 = V 



1+n 



1+n 2+n 
•2 2v 



V 



3 + n 2 + n~'~l + n 



In contrast to Cj(n, A, vq), 
form: 



Ci«(n,A) = \\^in + 4) + jE) 



the final expression for Cf{n, A) = C^(n, A, 0) takes a simple 

(C.12) 



18 - 40A - 12A3 + 2A'^ + SSA^ -6 + 4A 
n + 1 n 



-6 + 26A + ISA^ - 4A^ - SSA^ -SA^ + 4A - 2 + A^ 
-|- 2 A 



2 + n 



3 + n 



+A^ 



1 114A - lOSA^ + 36A=^ - 6A^ - 45 ^ 1 -86A + 129 



6 n 



+ 2 



2 n+1 
1 -66A + SOA^ - 44A3 + 12A'^ + 21 1 2A - ISA^ + 24A^ - 18A^ + 6 

^6 ~ 



2 + n 



3 + n 



on 



It follows, in particular, that the resummed width, integrated with respect to the lept 
energy (no cut, i.e. xq = 0), is: 

1 dr(A,r) _ f dn dr„(A) 

^ " {^.16) 



with 
dVniX) 



d\ 



To dXdr 

dTniX,x 



dX 



dx 



matched ■''1— A 



dXdx 



matched 



(C.14) 
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^g(A) + «(A) +C?(n, A)) ^^^iM + o(a^) 



TT 



X 



exp<^ 



Jo u 

-Bjiu)r{-u 



r(n) 



55(u)r(-2«) 
1 



Tin) 



r(n - 2n) r(l - 2u) 



T{n-u) T{l-u)J_ 
where Wq i{X) are given in Eqs. (C.2) and (C.4). 

D. Phase-space integrals for partially integrated width 
D.l Partially integrated width in hadronic variables 

We have seen that the resummed moments dTn{X,x)/dXdx of Eq. (3.40) (or Eq. (3.45)) 
where n is the Mellin conjugate to r = (P+ — A)/ (P~ — A) can be used to compute the 
differential width using Eq. (4.4). Similarly, the partially integrated width with a cut 
P"*" < P^ax can be computed by 



1 dr(p+ < p+,„ P-, El 

To dP- dEi 



mt 



L 



c+ioo 



dn 



1 - 



2m{n - 1) 

max 



P- 





\ n—l 




) ^ 


dTniX = 


P--A x-^^^) 
mi, ^ mi, > 


dXdx 



(D.l) 



This relation is useful in computing partially integrated widths with a variety of experimen- 
tal cuts by performing integration over Ei and P~ as well as the inverse Mellin transform 
and the Principal Values Borel integration numerically. 

The expression for the width with a cut on the invariant mass, in addition to a cut on 
the lepton energy, takes the form: 



V{P+P- <Ml,Ei> Eq) 



I 



Mb/2 i-Mb 

dEi e{Ei > Eo) / dP- 
JMb-2Ei 



(D.2) 



e{MB - 2Ei < M'x/P-) 



+0{M\/P- <Mb-2Ei) 



1 dr{P+ <Mb- 2Eu P-, El) 



dP- dEi 
1 dr(P+ < Mj/P-, P-, El) 
To dP~ dEi 



so Eq. (D.l) is used twice, in the latter term P^ax ^^t by the external invariant mass 
cut value M^/P~, while in the former by the actual phase-space limit Mb — 2Ei. Note 
however that the integrated distribution over the entire P+ phase space (P"*" < Mb — 2Ei) 
simply corresponds to the first moment. 



1 dr{P+ < Mb- 2Ei, P-, El) _ 2 dTn=i{X 
dP- dEi ~ ^ 



P--A 

mb 



mb ' 



dX dx 



(D.3) 



and therefore need not be computed by an inverse MclHn integral; such is necessary of 
course for the term where the upper limit on P^ is set by M^/P~. 
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When computing the phase-space integrals in Eq. (D.2) one should distinguish between 
three cases: for cuts where > — 2EqMb (in Fig. 4 this would correspond to 
< Eq) the hadronic invariant-mass cut is irrelevant and 

r{P+p- <Mj„Ei>Eo) /-^^/^^ /-^^ 1 dr{P+ <Mb- 2Ei, P-, El) 

= / dEi dP .(D.4) 

J-o Jeo jmb-2Ei io "-T ahii 

For cuts with {Mb - 2Eo)^ < M\ < M| - 2EqMb (in Fig. 4 this would correspond to 
Ef^ <Eo< Ef^y. 



T{P+P- <Ml,Ei> Eq) rh{MB-Ml/MB) 
= / dEi 



JEo 



MI/{Mb-2Ei) 

P- 

Mb-2Ei To dP dEi 



1/(^B-2E,)^^_ 1 ^r(P+ < Mb - 2Ei. P-, El) 



^ /-^^ 1 dT{P+ <Mj,/p-,P-,Ei) \ 

Jm%/(Mb-2Ei) To dP- dEi J 

+ r'" iE,r iP-L'^nP-<Ms-2E„p-.E,) ^^^^ 

J\{Mb-m\/Mb) Jmb-2Ei To dP dEi 

Finally, for cuts M|. < {Mb - 2Eq)'^ as in Fig. 4 where Eq < Ef'' < Ef^: 

T{P+P- <MI,Ei>Eq) _ rhi^B-M^)^^ .Mb 1 ^r(P+ < Ml/P-, P-, El) 

To Jeo Jmb-2Ei To dP- dEi 

^ .\{MB-MyMB)c rMl/(MB-2E,)^^_ ^ ^^p+ ^ Mb - 2Ei, P- , Ej) 

J^{Mb-Mx) I Jmb-2Ei To dP~ dEi 

^ /-^^ 1 dT{P+ <MI/P-,P-,EI) \ 

Jm'^/{Mb-2Ei) To dP- dEi J 

f^B/2 1-Mb I dT{P+ <Mb- 2Ei, P-, El) , , 

+ / dEi \ dP- ^ ,^ (D.6) 

J\{Mb-m\imb) Jmb-2Ei To dP dEi 



The last expression, corresponding to a stringent (charm-discriminating) hadronic-mass 
cut and a mild lepton-energy cut, is most useful for extracting \Vuh\- 

D.2 Distributions and partial widths with lepton energy cut 

The analytic integration over the lepton energy performed in Sec. 3.5 and Appendix C can 

be used as an alternative to the numerical integration over Ei . Starting with the rcsummed 
distribution of Eq. (3.52) with x > xq and converting to hadronic lightcone momenta using 
Eq. (4.1) we obtain the following expression for the double differential distribution: 

1 dT{p-,P+,Ei >Eq) _ 1 f dn _ KX~' (i -^y\ (D.7) 



dP-dP+ P- Jr-ioo 27ri V P~ J \ P~ 



P~-A ^ ^ ^„ _ 2Eo 

rrih dX 



1 drn{x=^,x>XQ='^) 
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The corresponding integrated width with a cut on is: 

1 driP-,p+<p+^,Ei>Eo) ^ r+'^ dn ( Ay-y. ^+axV:! 

7,_,^ 27r^(n-l) V P-) V P- ) 



dP- 



rrib dX 



(D.8) 



In order to compute the integrated width with an experimental PjJax cut one has to sum 
over two contributions, one from the small P~ region, where the experimental cut as well 
as the resummation are irrelevant, and one over the more important region of larger P~ 
where they are relevant. The result is: 

r{p+<p+^,Ei>Eo) ^ /-^-iax ^r„=i(A = ^ > ^) 

To Jo nibdX 

f 

JPr 



""^^ 1 dT{P-,P+ <P+^,Ei>Eo) 
lp+ To dP- 

max 



+ / dP- 



where in the second term Eq. (D.8) is used. Recall that the matching term in '^^^^^^^^o) 
has different analytic expressions depending on the sign of Mb — 2Eq — P~, see Eq. (3.53) 
and below. 

The double differential width of Eq. (D.7) can be readily converted into other kinematic 
variables, for example, the hadronic and leptonic invariant-mass distribution with lepton 
energy cut is: 



1 dr{M]„q^Ei>Eo) _ / ,p+,p_ r+-dn (^-7^)" / P+V" 
ro dM],dq^ J 27ri P' ^ P- ) 



X 



^rn(A = X > ^) 



5{Ml - P+P-) 5(g2 - (Mb - P+)(Mb - P")) 



mi, dX 

. n-l 



Lioo 2^ {P-?-Mj, [ {P^J 



dr„(A = ^,x>^) 



Mb mb dX 



Ml 



(D.IO) 



P- = Ej + ^fEf= 
E^ = {Ml-Ml-q^)/{2MB) 



Eq. (D.8) can also be used to obtain the partial width with a cut on the hadronic 
invariant mass M|-, i.e. 

V{P-P+<MI,E,>E0) ^ f^- dr.=i(A = ^,x>gg) 

To Jo rubdX ^ ■ ' 

+ dP- ^ ^^(-^"'-P^ < Ml/p-,Ei > Eq) 

Jmx ^0 dP~ 

where in the second term Eq. (D.8) is used. 
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